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Chapter 1 
Introduction 



1.1 Mobile ad hoc networks 

Mobile ad hoc networks (MANETs for short) are wireless networks which organize them- 
selves. They consist of a number of independent computing devices with wireless trans- 
ceivers. Not relying on a pre-existing infrastructure (such points or base stations), 
these devices can exchange data with each other in a multi-hop fashion, where each of the 
nodes is able to act as a router. A typical property of these system is that the position of 
network nodes, and hence the network topology, is not predefined, but is set up at random 
(or cannot be controlled) and will even change dynamically if the nodes are mobile. 

Such networks have been the subject of research since the beginning 1980s (usually 
under the term packet radio networks). But only recently, with the advent of commonly 
available, inexpensive wireless devices, the subject has gained much attention and focus 
[HHL03] . 

While MANETs are not currently in widespread use, there are a number of promising 
applications: Traffic information might be transmitted to cars on a motorway via an ad hoc 
network composed of on board communication systems [HBE + 01J. Pedestrians carrying 
mobile phones or PDAs might use a MANET for mobile data access, and firefighters could 
use such networks to gather critical information on scene |.TCH+04| . MANETs could prove 
to be particularly useful in situations where no network infrastructure is available, such 
as in disaster areas. Another interesting use case is their ability to extend an existing 
infrastructure, such as WLAN hotspots, beyond the range of the installed access points. 

A special case of MANETs, which is specifically the target of current research, are 
sensor networks ASSC02J which are formed of small, inexpensive, autonomous devices 
that are able to capture measurement data of some kind and transport them to a central 
location using ad hoc network techniques. Such sensor nodes are usually not mobile, but 
might be deployed in an area at random, e.g. by dropping them from a plane. Sensor 
networks might deliver valuable data for use in environmental monitoring, agriculture, or 
forest fire detection, just to name some examples. 

On the technical side, MANETs would typically be based on existing wireless technol- 
ogy, most notably IEEE 802.11 WLAN jXSOl] and Bluetooth |KRSW03j . While MANETs 
thus inherit a number of difficulties that are inevitably connected to wireless communica- 
tion - such as low reliability of links and the hidden-station problem -, there are a number 
of characteristics and issues that are specific to the ad hoc domain: 
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Due to the lack of a central controller, virtually all algorithms used on layer 3 and 
beyond must be distributed to the network nodes in order to provide for scalability and 
fault-tolerance. This applies not only to the application layer, but in particular to the 
routing protocols used. These are especially important since the network topology may 
change frequently at run time, so that traditional routing approaches cannot be applied as 
usual. A number of specialized routing protocols have been developed that are optimized 
for the specific needs of MANETs |Raj02| . 

Further, when using mobile or sensor devices, one is concerned with the problem of 
energy consumption: The network nodes are usually powered by batteries whose capacity 
may restrict the online time of each node (since users can only recharge their devices at 
intervals) or, in the case of maintenance-free sensor devices, may even limit the lifetime of 
nodes. A good part of the power consumption of the nodes is in fact due to their radio 
transceiver. Usual power-saving strategies are based on switching the devices into some 
inactive mode when not used; this is not favourable in ad hoc networks, however, since 
each device may be needed as a router. Thus, the range assignment problem is of central 
importance for MANETs: If the transmitting power (hence the radio range) of the network 
nodes is chosen too large, this amounts to a waste of battery resource. (It might also put 
limits to spatial channel reuse.) On the other hand, choosing the radio range too small may 
impact the network quality, since the number of point-to-point links is reduced. Therefore, 
a critical point in MANET design is to select the right radio range for a given density of 
nodes (or vice versa). Also, methods have been investigated to choose the radio range of 
each node dynamically jKKKPOOl iRRHOOj . 

Last but not least, a large number of research activities focus on the development 
of applications using MANETs, on specialized middleware (e.g. |Rot02| lHer03| ). and on 
security aspects BH03j . 

In evaluating design proposals for MANET systems, it is usually very hard to actually 
verify them in real measurements: Such experiments need to rely on a prototype implemen- 
tation, which is usually available only very late in the development cycle; moreover, they are 
quite cost-intensive, considering the large number of network nodes involved. Due to these 
obstacles, only quite few experimental evaluations have been performed MB.I00l lKNG+04j . 
with the MANET size being far below 100 nodes - which is rather on the low end for 
possible applications, considering that sensor networks of several 10.000 nodes are being 
discussed. Evaluation of protocols, etc. are therefore often based on numerical simulations: 
Using a statistical model that describes the spatial distribution of nodes and, for models 
with mobility, their movement on the deployment region, it is possible to evaluate the per- 
formance of routing algorithms, or to analyse general quantitative properties of MANETs; 
several off-the-shelf network simulators are available. 1 On the statistical side, a variety of 
different mathematical models are used to describe the mobility of nodes [CBD02J, one 
of the best-known being the random waypoint model. Recently, a number of publications 
have questioned the accuracy of the results of such simulations, both with respect to the 
properties of the statistical model involved |CSS02|. IYLN0 3J and to over-idealized mod- 
elling assumptions |KNG + 04] ; the accuracy of simulations must therefore be regarded as 
an open issue. Exactly solvable models, or other analytical results in MANET models, 
are currently only very rare, largely due to the complexity of the problems involved (see, 
however, the next section). 

1 Examples include OPNET (http://www.opnet.com/products/modeler/home.html), NS-2 
://www. isi.edu/nsnam/ns/), and GloMoSim (http://pcl.cs.ucla.edu/projects/glomosim/). 
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1.2 Quality and connectedness 

Let us take a closer look at the results known in the literature for measuring the quality of 
MANETs. Here we do not refer to the performance of routing algorithms or higher-level 
protocols, but we are rather interested in restrictions on the lower layers, related e.g. to 
inter-node connectivity. 

One natural question in this context is whether the MANET is connected, 2 i.e. whether 
there is a multi-hop network path between each pair of nodes in the MANET. Since we are 
dealing with nodes which are distributed at random, we are thus asking for the probability 
that the network is connected. 

Some early works |PPT89[IPir91j established asymptotic estimates for the probability of 
connectedness in 1-dimensional systems and conjectured analogue results for 2-dimensional 
systems. Here the nodes were distributed on an area (or line segment) according to a 
Poisson process of homogeneous density. The authors dealt with the probability that the 
area is completely covered by the MANET, i.e. that each point is in the range of at least 
one network node. Recently, the results for 2-dimensional systems were made precise by 
Xue and Kumar |XK04j . Denoting the number of network nodes by n, these results show 
that the mean local density of network nodes must to grow by a factor 0(lnn) in the limit 
n — > oo if one wants to keep connectedness (or coverage of the area). 

The probability of connectedness was also considered by Santi and Blough SB03J, based 



on earlier similar work [SB02, SBV01J. The authors derive asymptotic estimates mainly 
for the 1-dimensional system and present numerical (simulation) results also for 2- and 
3-dimensional systems. 

Bettstetter |Bet02j generalized this to the even stronger condition of /c-connectedness. 
(The network is called k-connected if between each node pair, there are at least k in- 
dependent network paths.) Using results from the theory of random graphs |Pen99j . he 
established analytical estimates for the 2-dimensional case and verified them with numeri- 
cal results. He also calculated the probability that none of the nodes is completely isolated 
in the network. 

More general quality measures have been defined by Roth |Rot03j and investigated in 
a numerical simulation. Here, not the probability of connectedness is taken as a qual- 
ity indicator (since, as the author notes, connectedness is a rather strong condition for 
MANETs); instead, measures based on the number of separated network segments, the 
size of these segments, and the dependence of these on changes in the network (e.g. a node 
being switched off), are being considered. 

Further, an analytical estimate of the bandwidth available to each node has been given 
by Gupta and Kumar [GKOOJ. The authors show that this bandwidth is of the order 
WJ y/n logn, where W is the bandwidth of the point-to-point links; thus, the throughput 
that each node is able to use rapidly decreases with the network size. 



1.3 Scope of this work 

In this work, we shall be concerned with the evaluation of quality measures for ad hoc 
networks described by statistical models. As in the last section, we will not consider 



2 Sometimes, the term strongly connected is used to describe this situation. 
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complex MANET protocols, but rather focus on simple models for connectivity between 
the network nodes; we shall establish explicit analytical results for the expected quality of 
MANETs on this level. 

We will set out from a statistical model of MANETs, similar to that considered in 
[SB03 , where a number of nodes is distributed independently at random in a given area. 
Focusing on the 1-dimensional situation (which might be interpreted as a network of cars 
on a road, or pedestrians on the sidewalk), we will show that the model is exactly solvable, 
and derive precise results for the probability of connectedness and other quality measures. 
Comparing these results to the literature, we will find that the numerical results both by 
Santi et al. Sl TVOTl and by Roth |E£t03j can be explained by our calculations, although 
they were based on different modelling assumptions. 

The work is organized as follows: 

In Chap. El we will define the general framework of statistical modelling that our 
calculations are based on. Idealizations and assumptions involved in this modelling will be 
discussed. 

Chapter El focuses on a specific model, the 1-dimensional MANET with homogeneously 
distributed nodes. Neglecting boundary effects, we will explicitly calculate the probability 
of connectedness for a fixed MANET size, and establish an asymptotic formula for the limit 
of large MANETs. This allows us to compare our results to existing work, in particular 
|SB03] . 

More general quality measures will be defined and analysed in Chap. EJ We classify 
these quality measures according to their scaling behaviour. Using the same model as in 
Chap. 01 we are able to establish explicit values for these measures in the 1-dimensional 
case, both at fixed size and in the limit of large systems. We compare our results to the 
numerical data from the literature [Rot03j and discuss similarities and differences. 

Chapter El discusses extensions of the results established in the previous chapters to 
more general situations. As an example, we explicitly treat a model where network nodes 
are switched off at random. An outlook is given to results in higher-dimensional systems 
and other extensions of the current results. 

Two appendices cover matters that are somewhat outside the main line of argument: 
Appendix 1X1 develops some mathematical results used in the main text, while Appendix iBl 
discusses certain issues found in the comparison of our results with SB03j. The reader is 
also referred to the index of notation on page 



Chapter 2 

Statistical Models for Ad Hoc 
Networks 



It is a characteristic property of ad hoc networks that, unlike in traditional infrastructure- 
based networks, the positions of the network nodes cannot be controlled. Therefore, it 
is generally useful to assume that the nodes are distributed at random in some area - in 
particular when the number of nodes is large -, and to use methods from probability theory 
in order to analyse the behaviour of the system. 

Such an analysis can, quite generally, be divided into two steps: 

(i) the definition of a mathematical model that represents the situation under discussion, 

(ii) the evaluation of this model and an analysis of its predictions. 

Both steps will, in general, involve simplifications and approximations of the "exact" 
situation: In the definition of the model, one decides on which aspects of the real situation 
should be modelled and which should be omitted; in the evaluation of the model, one often 
uses approximations (such as asymptotic expansions or limits) that only provide a certain 
level of precision. 

Certainly, the two steps are not independent: When defining the mathematical setup, 
one naturally has to take care not to define the model too detailed, in order to keep the 
complexity of evaluation within reasonable limits. So there usually is a trade-off between 
precision in modelling and precision in evaluation. 

However, it seems important to keep the distinction between the two steps clear, and to 
define clear interfaces between them. This is particularly important for the comparison be- 
tween different models or numerical approximations. It has recently been exposed (CSS02 
that the different simulation approaches can lead to very different results in the evaluation 
of protocols, even with regard to qualitative predictions. In this kind of situation, it would 
certainly be helpful to have a clear distinction between modelling and evaluation, since 
this might serve to clarify differences between the approaches, and to determine whether 
the difference lies in numerical approximation or in the general assumptions. It is also 
possible that an "informally" defined mathematical model (that is defined only by spec- 
ifying a numerical approximation) might include implicit properties that were not meant 
to be included in that way, as has recently been discovered with regard to the spatial node 
distribution in the random waypoint model YLN03J. 
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In this chapter, we will describe a general framework for the statistical description of 
MANETs, and discuss the assumptions and simplifications associated with it. We will try 
to define this framework quite generally, although only a very specific case will be analysed 
in detail in later chapters. This done is to provide a broader discussion of modelling 
assumptions and to hint at extension options for more complex systems. 

The formalism that describes the statistical behaviour of nodes is presented in Sec. 12.11 
while general assumptions related to the network model are discussed in Sec. 12.21 The 
evaluation of specific models is part of Chapters El to El 

2.1 A general statistical model 

We analyse an ad hoc network of n network nodes. These nodes are, at fixed time, dis- 
tributed at random over some volume or area. 

First, we will define the statistical side of the model, i.e. define the random location of 
nodes and, optionally, additional inner parameters. We will assume a sample space of the 
following form: 



Here fi spat iai denotes the sample space which describes the location of a single network 
node. It would usually be a subset of M. d , where d G {1,2,3}. The sample space fii nte rnai 
describes internal parameters of the node; e.g. the node might be switched off at a certain 
probability, or its transmission range might vary according to a random process. We are 
describing each of the n network nodes with the same sample space; note that this does 
not yet imply that the corresponding probability distribution is equal for each node, or 
independent between nodes. 

Another part of the sample space might be used to describe global random features 
of the model, e.g. the position of a shielding wall that disturbs network transmission. 
However, we will not make use of such an alternative here. 

On Q n , we then need a probability measure fi n which describes the distribution of 
nodes. We will specify further assumptions on fi n below. 

Let us now discuss the general modelling assumptions that are already implicitly in- 
cluded in the definition (|2.1jl . and additional assumptions that are typically made in order 
to simplify the discussion. 

Fixed number of nodes. With the above setup, we have assumed that the number of 
nodes in the network is fixed, i.e. it does not vary at random. This is certainly a restriction, 
since in a real scenario, the number of nodes might not be determined a priori; e.g. users 
might enter or leave the range of the network, or switch their devices off in order to save 
power. On the other hand, while it is possible to include a varying number of nodes into 
the setup of the sample space, this would increase the mathematical complexity of the 
system considerably, since it would require to move from a usually finite-dimensional space 
or manifold fl n to an infinite-dimensional situation. 

In our context, it seems to be justified to stay with the situation of fixed n for two 
reasons: First, when considering a large number n of nodes, one expects that the effect 
of a varying number of nodes is small, and that it suffices to take only the mean number 
into account. Second, if we explicitly need to account for nodes dynamically joining the 
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network, we can always model them as nodes which are randomly switched off from the 
network, including this aspect as a feature of fi interna i. Such an analysis will be presented 
in Section EIU 

Note that while we model the statistical situation only for fixed n, we will usually be 
interested in the expectation values of random variables "for large n," i.e. in the limit 
n — > oo. 

Static situation. Our analysis restricts to the situation of the network at fixed time. 
This may seem to be a bit contrary to our goal to describe mobile ad hoc networks. 
However, mobility of the nodes does not imply that the probability to find a node within 
a specific region varies with time. In fact, since the movement of nodes (e.g. of visitors 
in a shopping centre) will usually not be under our control in realistic situations, the best 
assumption might be that at any fixed time, nodes are distributed according to a static 
probability distribution. 

In other approaches to the statistical description of MANETs, one often introduces an 
explicit model for the random movement of nodes (such as the random waypoint model). 
However, even in these models, one would assume that the spatial distribution of nodes 
stays constants over time, or rather consider it as a problem in the model if this is not the 
case |YLN 03j. In fact, network simulators may need a "warm-up phase" until a "steady 
state" in spatial distribution is reached. 

Regarding time averages of random variables, we can deduce results from our static 
model if the network system is ergodic: This means that time averages can be replaced 
with averages over the spatial coordinates of nodes at fixed time, which we can handle 
directly. Equivalently, we may require that for each initial configuration of the network, 
we reach almost every other possible configuration after waiting for a sufficiently long time 
(Birkhoff's ergodic theorem; see |Pet83j). Ergodicity of the system is not guaranteed and 
depends on a mobility model still to be chosen; however, lacking specific information on 
the time dependence of the system, it seems to be a natural assumption for our purposes. 

Certainly, our model could easily be extended to describe non-static situations by mak- 
ing the probability measure dependent on the time t. However, our setup generally does 
not allow to describe aspects of the system that involve direct time dependence of random 
variables. For example, assuming ergodicity, we might be able to answer the question: "For 
what portion of the time is a specific node connected to the network?" , but our setup does 
not allow to discuss the question: "For how long does a specific node stay connected to 
the network, once it has established a connection?" In our discussion of quality measures, 
such time dependencies will not be relevant; for a discussion of routing algorithms, on the 
other hand, they may be a crucial feature. 

Independence. In addition to the above assumptions, we will apply another simplifi- 
cation, namely the statistical independence of the nodes. On the mathematical side, this 
means that our probability measure is reduced to a product 



n 



(2.2) 
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where /ij are probability measures on f2 spa tiai x ^internal, describing the distribution of a 
single node. 

On the modelling side, this means that the different nodes will have no mutual influence 
on their positions (or other internal states). This need not be fulfilled in realistic situations; 
for example, in a traffic jam, the position of a specific car will very well be influenced by 
the position of the car in front of it. Such aspects cannot be described when making the 
above assumption; however, it seems plausible that in many situations, such effects will 
not play a major role. 

Identical distribution. In addition to the independence of nodes, we will assume in all 
our examples that the node are distributed identically, i.e. that all /z° ode are in fact equal: 

Vj : /i° odc = /i^ de (2.3) 

This seems to be a natural assumption if there is only one type of node involved in the 
network. It does not cover a situation where certain nodes are distinguished from others, 
e.g. where certain users prefer a specific part of the area. We might, however, still cover 
these situations when modelling this behaviour within ilmtemai; that 1S , we would let the 
user choose "at random" which area he prefers, while preserving the identical distribution. 
(However, such models will not be covered in this text.) 

No feedback. In all of the following text, we will assume that the probability measure 
H n is given a priori as a fixed quantity, and that it does not depend on the details of the 
MANET quality; alternatively speaking, there is no "feedback" from the random variables 
to the probability distribution. 

Of course, one might in principle think of a situation in which users prefer to visit 
areas where the MANET quality is usually good, or in which they tend to switch off their 
devices if they loose connectivity for an extended period. Such aspects would need to be 
modelled in form of a (supposedly complicated) relation in \i w e.g. a differential equation, 
that would leave us with the task of finding a solution for ii n before calculating expectation 
values. However, this lies far beyond the scope of the current presentation. 

2.2 Random variables 

Having specified the statistical behaviour of the system, we will now turn to a description 
of the random variables. Random variables would include, e.g., the number of network 
segments, the number of nodes that a specific node is connected to, or the spatial dis- 
tance between two nodes. In the general mathematical setting, a random variable is an 
(integrable) function 

F : a n -> R. (2.4) 
As usual, we consider the expectation value of F, defined as 

E[F] := I dfi n (u)F(u). (2.5) 
and interpreted as the statistical mean of F. We sometimes write it as F for short. 
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Without putting too much emphasis on the mathematical formalism, it should be noted 
that a random variable itself does not include the statistical description of the model; for 
each fixed u G fi, its value F(u) is simply the "deterministic" value of the function in the 
elementary event u. The statistical behaviour is described via the expectation value alone. 

Usually, the definition of random variables will depend on n, as well as on other param- 
eters of the system (such as the range r of the radio devices). As above, we will often not 
denote this dependence explicitly, in order not to overburden the notation; in case were it 
becomes necessary, we will explicitly write F^ n \ F^ n,r \ or similar. 

In order to fix our notation, let us briefly introduce some special random variables, 
which are connected to events on Q n . An event EV is a subset of the sample space fl n ; 
as an example, take the event CONN which contains all points of Q n that correspond to 
situations where the MANET is strongly connected. For notational purposes, we will often 
write events as Mev C fi„ when referring to it as a set. To each such event corresponds 
its characteristic function Xev, defined as 



is the probability that the event EV will occur. 

After these formalities, let us discuss our modelling assumptions on the random vari- 
able side more closely. It is difficult however to investigate properties of specific random 
variables (such as connectivity of the network) without specifying a concrete model, which 
we postpone to Chap. H3 However, we shall discuss a number of general assumptions on 
the random variables, and how we wish to handle them. This follows a recent discussion by 
Kotz et al. |KNG + 04] who identified a number of common assumptions in MANET models 
and compared them with experimental results. The authors criticized these assumptions as 
begin too restrictive for realistic scenarios; we will in fact stick to all of these assumptions 
in this text, and will argue in the following why they are justified in our simple situation. 

The world is flat. While radio propagation is a 3-dimensional phenomenon, the nodes 
of a MANET are usually distributed over some 2-dimensional (e.g. sensors deployed in an 
area) or even 1-dimensional region (e.g. cars on a road). Truly 3-dimensional situations 
will only very seldomly be found in practice, since ceilings in buildings, etc. usually block 
radio propagation. If some network nodes are located in vertically exposed positions (e.g. 
on hills), it seems more appropriate to include this effect in the model by modifying their 
individual radio range (see below) rather than turning to a 3-dimensional description of 
their position. In most of this work, we shall restrict to the 1-dimensional case for simplicity. 

A radio's transmission area is circular. 2-dimensional MANET models usually as- 
sume that the range of network nodes is not dependent on direction. While this seems to be 
a natural assumption at first, it is often not realized in experiment |KNO + oH Fig. 1]; in par- 
ticular, commonly used antennas are not omnidirectional. However, for the 1-dimensional 
situation that we will consider, these properties will obviously be of less importance. 




(2.6) 



which is a random variable in our sense 



Its expectation value 



Pev := E[xev] 



(2.7) 
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Signal strength is a simple function of distance. In generalization of the last point, 
it may even be difficult in experiment to find any simple relation between the spatial 
distance of nodes and the signal quality on point-to-point links, since the signal strength 
is influenced by radio reflection, shielding obstacles (including e.g. the person carrying a 
mobile device), and other effects that are difficult to control. In fact, the data presented in 
|KNG + 04j suggests that the radio range of nodes should rather be described by a statistical 
process. We might include this behaviour in our model by assigning the radio range of 
nodes at random, albeit at the cost of a much increased complexity in evaluation. However, 
for the simple connectivity properties we will consider, it seems reasonable that only the 
mean radio range of nodes will be relevant for our results - see also the discussion in 
Sec.O 

All radios have equal range. We will assume in our specific models that all nodes are 
equal with respect to their radio range. Due to varying background noise, differences in de- 
vice configuration, and also for reasons named above, this may not be given in experimental 
situations. Again, we might include this in our model by considering the radio range of 
nodes as a random variable, or choosing it dependent on the node's spatial position; we 
will however refrain from doing so in the present work. 

If I can hear you, you can hear me (symmetry). Many MANET protocols discussed 
in the literature rely on network links to be bidirectional, while it has been stated in 
|KNG + 04j that this assumption is often not valid in practice; in particular, packet collisions 
may lead to unidirectional links. While this may be a crucial feature for routing protocols, 
unidirectional links should not affect our simple evaluations of connectivity: We aim at a 
description of the connectedness between nodes and disregard packet loss rates, etc. 

If I can hear you at all, I can hear you perfectly. For our evaluation of MANET 
connectivity, we will focus on the question whether a point-to-point link between two nodes 
can be established, and will not aim at a calculation of network throughput, packet loss, or 
error rates. Therefore, we can assign a sharply defined "range" to each node, below which 
we assume point-to-point links to be established, and beyond which no communication 
is possible. Certainly, for a more detailed analysis of MANETs, it should be taken into 
account that there is no sharp spatial cutoff for connectivity, but that the signal strength 
decays gradually with increased distance from a node. 



Chapter 3 

The Connectivity of 1-dimensional 
Networks 



We will now proceed to a specific MANET model which we will analyse in detail. This 
model restricts to the 1-dimensional situation, i.e. the nodes are deployed at random along 
a straight line. One might think here of pedestrians moving along sidewalks, or of cars on 
a road, that carry wireless devices. While this model is relevant at least for parts of the 
proposed applications, it turns out to be particularly simple in mathematical description, 
so that we can derive explicit analytical results e.g. for the probability of connectedness. 

In Sec. 13. 11 we will first give a definition of the model and introduce specific assumptions. 
Then, in Sec. 13.21 we consider a variant of the model - the model with periodic boundary 
conditions - which allows us to calculate the probability of connectedness, both for a 
fixed node number n and in the limit n — > oo. Section EPl will then return to the model 
with "usual" boundary conditions, transfer our results to that situation, and compare the 
outcome with analytical and numerical results known in the literature. 

3.1 Definition of the model 

As mentioned in the introduction, we will now consider a 1-dimensional system with n 
network nodes distributed at random. More specifically, we assume that the network 
nodes are distributed on an interval [0, £], that is, we set 



where the space ^internal is trivial, i.e. we consider no additional internal random parameters 
of the network nodes. For simplicity, we will assume that the n nodes are distributed 
identically and independently, according to the equal distribution on f2 spatial . This means 
that 



This fixes the statistical behaviour of the system. It still remains to define the random 
variables of interest. 

In this chapter, we will mainly be concerned with the question whether the network is 
connected, i.e. whether all nodes are able to communicate with each other in a multi-hop 



Spatial =[0,4 n n =[0,£] n , 



(3.1) 



(3.2) 



- 15 - 
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fashion. To that end, we will assume that all nodes have a fixed (and identical) radio range 
of r. Two nodes with coordinates X{ and Xj can communicate directly with each other if 

\xi — Xj\ < r. (3.3) 

It is then clear what "connectedness" means in the model. This is essentially the situa- 
tion considered by Santi and Blough SB03J, who derived estimates on the probability of 



connectedness in the limit £ — > oo. 

Before we proceed to the precise definition of random variables and the calculation 
of expectation values, let us first discuss some general properties of the random variables 
involved, since the model has some symmetries that we will exploit to ease our calculation 
later on. 

The first of these properties relates to the fact that the behaviour of the system does not 
depend on £ and r explicitly, but that is stays the same when r, £ and all coordinates are 
scaled by a common factor (there is no "fixed length scale" in the system). Heuristically, 
this is easily understood directly from the model; since we will make a lot of use of this 
property, let us however describe it more formally. 

Definition 3.1. In the 1- dimensional MANET model, a family of random variables 

F (n,i,r) . [Q^n _^ R ^ calM scaling % j 

Wn e N VA > Wx G [0,£} n : F (nAr \x) = F (n ' xe ' Xr) (Xx) . 

In fact, all random variables considered in the following will be scaling; this is due to 
the fact that the connectivity between two nodes is not affected by scaling, cf. Eq. ()3.3|) . 
The consequence of this property is that expectation values are indeed dependent on the 
ratio r ji only: 

Proposition 3.2. Let F^ n ' e ' r ^ be a scaling family of random variables. Then we have for 
allr>0,£> 0: 

Proof. By definition of the expectation value, we have 

Ej^nAr)] = J d n xr n F (nAr \x)=r n [ d n xF {n > 1 > r/e) (r 1 x), (3.4) 

J[0,£] n J[0,£] n 



using the scaling property with A = £ l . Now a simple substitution of variables x\= £ l x, 
leads us to 



E [ j p(n,<>r)] = f d n x' F^ l ' r/t \x') = E[F (n ' 1,r/<) ], (3.5) 

as proposed. □ 

Since in what follows, our results will relate to expectation values or probabilities, they 
will therefore depend on n and rji only. Alternatively speaking, we can refer to the sample 
space fl n = [0, l] n at any scale and use the "normalized radio range" p := r/£ in place 
of r, thus reducing the number of parameters by one. Since all our random variables will 
be scaling, it is justified to use this simplified model only; we will return to the explicit 
parameters £ and r only for comparison with experiment or other publications. 

The next general property is related to the fact that all n nodes are treated as equal in 
the model. 
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Definition 3.3. In the 1- dimensional MANET model, a random variable F is called sym- 
metric if, for any permutation a : {1, . . . , n} — ► {1, . . . , n}, it holds that 

Va? E tt n : F(xi, ...,x n )= F(x a{1) , . . .,x a{n) ). 

All random variables we consider in the following will be symmetric. 1 This property 
has the consequence that we can calculate expectation values more easily: In the integral 



E[F] = / d n xF(x), (3.6) 
J[o,i] n 

we can split the integration region [0, l] n into n! regions where the coordinate values are 
sorted in a specific order, i.e. R\ = {x | x\ < x 2 < ••• < x n }, R 2 = {x \ x 2 < x\ < 
£3 < . . . < x n } etc., ignoring sets of volume zero. Since all these regions have identical 
volume, and since a symmetric random variable F is not affected by a change in the order 
of variables, one has in this case 

E[F] = n\ I d n xF{x). (3.7) 



More explicitly, we can express this as 



E[F]=n\ f dx x I dx 2 ... ( dx n F(x); (3.8) 

Jo J x\ J x n —i 

this form is often convenient, whenever F can be formulated easier in the "sorted" coordi- 
nates X\ < . . . < x n . 



3.2 Connectivity with periodic boundary conditions 

Up to now, the system we defined was identical to the one considered by Santi and 
Blough [SB03]. In this model, one would define that in sorted coordinates x\ < . . . < x n , 
the node i is connected to its neighbour i + 1 if Xi + \ — Xi < p; for the nodes 1 and n, how- 
ever, there is no left-side or right-side neighbour, respectively, which they could connect to. 
While this definition seems somewhat natural, it leads to an increased complexity if one 
wants to derive analytical results: It includes a description of the effects at the boundary 
of the network, which one implicitly has to account for in any calculations. 

As a method to overcome these difficulties, we will introduce periodic boundary condi- 
tions in our model: We will say that the left-most node is connected to the right-most one 
if 

xi + 1 — x n < p. (3.9) 
fact, for any random variable F we might always define the symmetric random variable 

^Symm(a;) := — ^ F(x CT (i) , . . . , X CT („) ) 
a 

due to the symmetry of the underlying integration measure, we easily find ~E[F] = E[Fs ymm ]. 
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This amounts to a periodic extension of the node coordinates to the region outside [0, 1]. 
One might also think of the nodes being located on a closed path rather than an interval. 2 
While this change seems to be a bit technical, it is justified for two reasons: First, we 
are interested in the behaviour of the MANET in the "bulk" and not at the boundaries; 
it is thus reasonable to eliminate boundary effects via the periodic extension. (In fact, in 
a realistic scenario such as the shopping center example considered by Roth RotOS], the 
paths that users are located on would include both closed curves and open segments, and 
thus a "disconnected" boundary condition is a priori not more realistic than a periodic 
one.) Second, and more importantly, it is expected that in the limit of large MANETs 
(n — > oo), these boundary effects play no role, and both models lead to the same results. 
We will explicitly show this for the probability of connectedness in Sec. K13.ll 

3.2.1 Transformation of the probability space 

On the analytical side, the introduction of periodic boundary conditions amounts to a 
change in the random variables (the probability distribution is unchanged); it results in 
the following property. 

Definition 3.4. A random variable F : [0, l] n — > R is called translation invariant if 

Wx e [0, l] n VA G R : F(xx, ...,x n ) = F(x 1 + A, . . . , x n + A), 

where the function F is taken to be periodically continued to R n ; i.e. F(xi + l,x 2 , . . .) = 
F(xi, x 2 , ■ ■ .) etc. 

We will later see why all relevant variables in our context are in fact translation in- 
variant. Let us first analyse the consequences of this property. To that end, let F be a 
symmetric and translation invariant (as well as scaling) random variable. Its expectation 
value is given by Eq. (J3.8J) . In that integral, let us introduce the next-neighbour distances 
Hi = Xi+i — x-i (i — 1, • • • , n — 1) as variables; this results in 

/•l rl-Xi rl-xi-m /"i-Zl-ESTi 2 Vi 

E[F] = n\ / dxi / dyi I dy 2 ... dy n -iX 

Jo Jo Jo Jo 

n-1 

x F(xi, xi + y h x 1 + y 1 + y 2 , . . . , x x + } J Vi)- (3.10) 

i=i 

In the argument of F, we can certainly replace x\ with due to the translation invariance 
of F . Moreover, we set 

n-1 

F(y 1 ,...,y n ) = F(0,y 1 ,...,^2y i ), (3.11) 

where the purpose of the apparently "redundant" variable y n is as follows: If we set y n = 
1 — Xir^i 1 Vi-> then it is easily seen from the symmetry and translation invariance of F that 

2 The use of periodic boundary conditions is a well-known technique for dealing with similar types of 
boundary problems; it has also been applied the analysis of MANET connectivity before |Bet02| . 
3 Note that the definition does not refer to sorted coordinates. 
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F is shift-symmetric in the n variables, in the sense that 

F(y 1 ,...,y n ) = F(y 2 ,...,y n ,y 1 ). (3.12) 

Regarding the integration domain in Eq. ()3.1()|1 . we can see that the combined integration 
over x±,yi, . . . , y n -i runs over the n-dimensional standard simplex V n ; thus 

E[F]=n\ [ dx 1 d n ~ 1 y F{y). (3.13) 

JVn 

(The standard simplex and its properties are discussed in Appendix IA.1| which we will 
frequently refer to.) Choosing a different coordinatization of the simplex, we can express 
this as 

E[F}=n\ [ d n ^y [ ' = * * dx x F{y). (3.14) 

JVn-! JO 

The integration over x\ can easily be executed: 

/n— 1 
d n - l y (i-J2yi)Hy)- (3.i5) 
-n-l i= l 

Setting y n — 1 — Y^i=x Vh an d comparing with Eqs. ()A.12jl and ()A.13j) in Appendix IA.1| 
we can rewrite this as an integral over the top surface T n of the simplex in n dimensions: 

„ n 

E[F] = n\ / d n y5(l - V^) y n F(y). (3.16) 
J M n i= i 

Now noting that the integration measure is completely symmetric with respect to an ex- 
change of variables, and using the shift-symmetry of F [cf. Eq. f!3.12j) ]. it is clear that we 
can replace the factor y n in the integrand with any other yi without changing the integral's 
value; so we can as well replace it with the mean: 

„ n n 

E[F] = n\ / cfV(l-y>) -(T,Vi)F(y). (3.17) 
However, under the integral, we have Y17=iVi = ^- Thus, our result is 

„ n 

E[F] = (n - 1)! / d n yS(l - Ty,) F(y). (3.18) 

Comparing with Proposition IA.3| we can rewrite this as 

E[F] = [ dt£-*(y) F{y). (3.19) 

7[0,l] n 

the next-neighbour coordinates are distributed equally (not independently!) over the top 
surface T n of the standard simplex. Let us summarize: 
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Theorem 3.5. Let F be a symmetric and translation-invariant random variable on Q n = 
[0, l] n , considered with the equal distribution. Let F be the corresponding random variable 
[see Eq. 113. lip ] on Q! n = T n , considered with the equal distribution on T n . Then 

E[F] = E[F]. 

In fact, it will be more convenient in most cases to define the random variables di- 
rectly in terms of the next-neighbour coordinates; given that the so-defined variable F is 
shift-symmetric, we can always define an underlying symmetric and translation-invariant 
random variable F. We will not even distinguish the two associated random variables in 
notation (where this is unambiguous). 



3.2.2 Connectedness 

We will now turn to calculate the probability that the MANET is connected. This 
needs some explanation with regard to the periodic boundary conditions: We will call 
the MANET connected if all next neighbours are connected, including the left-most and 
the right-most one (which are connected "via the boundary"). More formally, we define 
the event CONN-PB in next-neighbour coordinates as 

Mconn-pb := {v e T n | Vj : Vj < p}. (3.20) 

We also consider the more general event k-DISCONN-PB for k e No, defined as 

M k _ DISC oNN-PB := {y e T n | yj > p for exactly k values of j}, (3.21) 

meaning that the network is disconnected at k places (or, equivalently speaking, into k 
segments). Note that CONN-PB = 0-DISCONN-PB. 

Our task is to calculate the probability of k-DISCONN-PB. A central tool for this is 
the inclusion-exclusion formula (see Appendix IA.2J) : it gives us 

Pk-mscoNN-PB = ]T(-ir * (fy Sj, (3.22) 
j=k \ ' 

where 

S,= P (y™ >pA...Ay w , >p). (3.23) 

{mi ,...,m,j} 

It remains to calculate the probability of the event under the sum, which is handled in the 
following lemma. 

Lemma 3.6. Let j e {0, . . . , n}, and let {mi, . . . , mj} C {1, . . . , n} be a j-element subset. 
Then 

\l-jpT- 1 ifj<l/p, 



P(y mi > P A . . . A y mj >p) 



otherwise. 



Proof. Let P be the probability in question. We will prove the result by induction on j. 
For j = 0, it is obvious that P = 1 as proposed. So assume that we have verified the result 
for j — 1 in place of j. The case j > 1/p is obvious, since Yli=i Vi = 1> so let j < 1/p in 
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the following. The characteristic function of the event can be expressed as a product of 9 
functions; 4 that results in 

P= I dnl~ c \y) f[e(y m -p). (3.24) 

J M n i= i 

Applying Lemma IA.4I with respect to the variable y mj , we obtain 



p=(i-pr- 1 [ d£-«(v) n% 

J m n i= i 



i-p' 

pT^PiVrm > ~r— A ... A y mj _ Y > -?—). (3.25) 
1 — p 1 — p 



Here we can apply the induction hypothesis for j — 1 in place of j and p/(l — p) in place of 
p; note that the condition j < 1/p guarantees that j — 1 < (1 — p)/p. This shows us that 

P = (1 - p) n - 1 (1 - (j - l)-^)^ 1 = (1 - jp) n -\ (3.26) 

1 — p 

which proves the lemma. □ 

Applying this lemma in Eq. ()3.23|) . and then inserting into Eq. ()3.22j) . we can establish 
an explicit expression for -Pic-disconn-pb- Note that in Eq. ()3.23|) . all summands are in fact 
equal, so that we only need to count the number of terms, which is (™) . Our result then is: 

Theorem 3.7. In the 1- dimensional MANET with periodic boundary conditions, one has 
for each k E No, 

[i/p] 



In particular, 



Pk-DISCONN-PB = ^(-l) J " fc Q (j) (1 " JP) 

j=k 

[1/P] XV 

PcoNN-PB=Y,(- i y( n Y i -jp) 



n-1 



n-1 



3=0 K " y 

Here [1/p] is the Gauss bracket of 1/p, i.e. the greatest integer which is less or equal 
to 1/p. Note that the formula is valid for [1/p] > n as well, since the factor (") evaluates 
to for j > n, so that these summands automatically vanish. 

We have thus found an explicit expression for -Pk-DiscoNN-PB; the expression is defined 
piecewise as a polynomial in p of degree n — 1. In particular for small values of p, the 
sum involves terms of high modulus and opposite sign; thus a numerical evaluation with 
floating-point techniques may lead to problems due to round-off errors. However, inserting 
p as a fraction, we can use integer arithmetics in order to evaluate the sum, thus bypassing 
the problems mentioned. 

Figure EIU shows the behaviour of -Pconn-pb, plotted against n (on a logarithmic scale) 
and np. Two things are noticeable: First, at fixed n, we obviously have -Pconn-pb - ► 1 
for p — > oo and Pconn-pb — ► 1 for p — * 0. This is expected and can directly be seen from 
the arithmetic expressions. Second, it seems that in the limit or large n, the probability 
-Fconn-pb is basically a function of one parameter np — In n. This asymptotic behaviour 
will be discussed in the next section. 



J See Eq. (| A. 3|> in Appendix IA. II for the definition of the Heaviside 9 function. 
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Figure 3.1: Probability of connectedness for the 1-dimensional MANET 
3.2.3 Asymptotic behaviour 

Apart from the probability of connectedness for fixed parameters p and n, we are partic- 
ularly interested in the behaviour of our model for large MANETs, that is, in the limit 
n — > oo. However, although Fig. 13.11 suggests that there is some well defined large-scale 
limit of the system, it is not apparent from Theorem 13 . 71 how -Pk-DiscoNN-PB behaves in this 
limit. In this section, we will discuss -Pk-DiscoNN-PB in the large-scale limit and derive an 
asymptotic approximation formula. 

Since the detailed calculation turns out to be quite technical, let us first present a 
heuristic sketch of the underlying ideas, where we will restrict ourselves to -Pconn-pb- We 
can rewrite the expression from Theorem 13.71 as 

[i/p] , 

Pconn-pb = ]T(-iy ' - {l-j p y-\ (3.27) 

Using Sterling's formula (Inn! w nlnn) and Taylor expansion (ln(l — x) ~ — x), we see 
that for large n and moderate j, 

ln_ T7 — ~v\ ~ i( lnn - ln -?)' ln ( 1 - ip) n 1 ~ -if™! ( 3 - 28 ) 

so the polynomial factor (1 — 2P) n X dominates the binomial factor for medium to large j, 
such that only a very limited number of summands (j < jo) will actually contribute to the 
sum in Eq. ()3.27j) . For these terms, we can individually let n — > oo at fixed j. Here we 
have 

— = n{n - 1) . . . (n - j + 1) w n> and (1 - jp) 11 ' 1 w e~ lpn . (3.29) 

(n-j)l 
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Inserting into Eq. ()3.27j) . this means that 



30 ' 



Pconn-pb « J2 —T^e- np y. (3.30) 

3=0 J ' 

Without changing the value of the sum significantly, we can replace jo with oo here; then 
the sum becomes an exponential series, and we see that 

-Pconn-pb ~ exp(ne~ np ) = exp(-exp(-np + Inn)). (3.31) 

Of course, controlling the limit n — > oo is in fact not as easy as suggested above, and 
we have to turn the heuristic arguments into a rigorous proof in order to be sure about the 
large-scale behaviour of the model. This is the content of the following theorem. 

Theorem 3.8. Let (p n ) be a sequence in R + , and suppose there is an r] E R such that 

, n— >oo 

np n — Inn > rj. 

Then, we have for every k G No' 



^k-DISCONN-PB > i\ eX Pl e 



-r\k 



->l\ 



Proof. First, let us note some properties of the specified limit: Since np n — Inn — > rj, we 
certainly have p n ~ \nn/n and thus 

p n -> 0, np n -> oo, npl -> 0. (3.32) 

Now let us turn to -Pk-DiscoNN-PB- We can rewrite the expression from Theorem 13.71 as 



n k e knpn , , ■ . n\n k 



k-DISCONN-PB 



fc! 7 (n-j)l(j-A;)! 

J— K 



Shifting the summation index by —k, this is equivalent to 

k -knp n i l IJ^ k n\ n - k 
Pk-DISCONN-PB = TT— Y\ ^7 — ^T7 (l " 0' + k)^ 1 ^ . (3.34) 

k\ z — ' ?! (n — 7 — /c ! 



In the factor that precedes the sum, it is obvious that 

n k e -knp n = e -fc(np n -lnn) e ~ V k . (3 35) 

so it remains only to control the convergence of the sum itself. We will next investigate 
how fast the summand terms a 3 - vanish for large j. We can certainly say that 

__ = n f n - 1) • • • (n - j - k + 1) < n 3+k (3.36) 

{n-j- k)\ 
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and thus 

lnaj < jinn + (n — 1) ln(l — (j + k)p n ) + /cnp ra . (3.37) 

Since it is known from the Taylor series of ln(l — x) that ln(l — x) < —x for all x G (— oo, 1), 
we see that 

lncij < j(lnn - np n + p n ) + kp n . (3.38) 
Now since np n — Inn — > n and p n — > 0, we can certainly find no such that 

Vn > n Vj : ma, < 2r?j + 1 or, equivalently, ctj < e V3+ . (3.39) 

According to Stirling's formula (cf. Theorem IA.7I in Appendix IA.3)1 . we can say that for 
any j 

>3_\j 

and thus for n > n 



j\ > (I)', (3.40) 



<e-(—y. (3.41) 
J ! 3 



Given e > 0, we can thus find j such that 

a 



Vn>n Vj>j : ^ < (3-42) 
3- ^ 



This means that 

[l/Pn]-fe 



Moreover, after possibly increasing jo, we can achieve that 

|^L_L e -,i|< 6) (3.44) 

since the exponential series Ylj^/j^- converges absolutely on R. Further, we can assume 
that 1/ p n > jo + k for n > n . 

It remains to estimate the convergence of the terms for j < j in Eq. ()3.34|) . To that 
end, note that the above estimates are uniform in n: Once we have fixed jo and no for given 
e, we can consider the limit n — > oo without changing jo. Thus, there are only finitely 
many terms left to estimate, and we can consider the limit in each of them individually: 
We want to show that for each j < j , one has 

a ./ e -nj i. (3.45) 

Explicitly, we know that 

a i/ e ~ Vj = r— >m ni C 1 " U + k)p n ) n - l e kn ^. (3.46) 
(n - j - fcj! 

Certainly, the first factor converges as n — > oo: 

raj n~ fc ~ J = n(n - 1) ■ ■ ■ (n - j - k + 1) n^oo ; , 4 s 

(n - j - fc)! ~~ * ^ ' ) 
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Furthermore, we see that 

In {n 1 (1 - (j + k) Pn ) n - l e knp ^) = j Inn + (n - 1) ln(l - (j + k)p n ) + knp n + nj. (3.48) 

Again, we use the Taylor expansion ln(l — x) — —x + 0(x 2 ); this results in 

dHHD = j(lnn - n Pn + n) + 0{p n ) + 0{np 2 n ). (3.49) 

According to Eq. (J3.32|) . all of the terms on the right-hand side vanish in the limit; this 
proves Eq. f)3.45j) . Since we had seen in Eq. (|3.39J) that the cij are uniformly bounded in n 
(at fixed j), we have a forteriori that 



P -m 



dj -e-" 3 \ < \aj\ |1 1 ^-^0. (3.50) 



This means that we can find rii > n such that for any n > n 

' 0_1 < -1)3 (-1)3 



i=o J i=o J 

Now combining Eqs. ()3.43|) . ()3.44|) . and ()3.51|) . we know that 

[l/ Pn ]-k . OO ■ 

Ve > 3n t Vn > n x : | ^ L -y^% - X] ^-^( e ~' 7 ) J | ^ 5e - ( 3 - 52 ) 

i=o ] ' j=o ] ' 

Rewriting the exponential series as an exponential function, this means 

[l/pn]-k , 

£ UL aj exp(-e^). (3.53) 

Inserted into Eq. ()3.34|1 . this proves the theorem. □ 

Let us add another result for the limit n —>■ oo, which we state for Pconn-pb only: 
Suppose that np n — Inn — > oo in the limit. Then for given rj, we can certainly construct a 
sequence (p' n ) with p' n < p n such that np' n — \nn n. Since -Pcqnn-pb * s mon °t° nous i n P 
at fixed n, we see that 

p(n,Pn) > p(n,p' n ) n-^oo / -7A /o r 4 \ 

-"CONN-PB — -'CONN-PB t^Pl e J- I^O.J^J 

We can choose arbitrarily high here; that means ^conn-pb ~~ ► 1- A similar result for 
np n — Inn — > — oo can be obtained in the same way. Let us note this for reference: 

Theorem 3.9. Let (p n ) be a sequence in IR + , and suppose that 

np n — Inn ► +oo or np n — Inn > — oo. 

Then, we have 

p conn-pb "^°°> 1 or > respectively, Pconn-pb > °- 
(A similar result could be proved for Pk-DiscoNN-PB, but we will make no use of it.) 
Figure 13.21 shows the quality of the asymptotic approximation of -Pconn-db for growing 
n. As might be expected from the details of the proof, the convergence is particularly fast 
for large n. For n < and low n, the absolute error is small, but on a relative scale, the 
approximation is rather unusable. This may be understood from the fact that Pconn-db is 
exactly for p < 1/n, as is apparent from the model, while the approximation exp(— e~ v ) 
still gives positive, if very small, values. 
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AP = exp(-e^) - Pconn-pb Q = p conn-pb / expf-e^) 




(a) absolute (b) relative 

Figure 3.2: Quality of the asymptotic approximation for Pconn-pb- Relative and absolute 
comparisons are plotted against n (on a logarithmic scale) and rj = np — Inn. 

3.3 Connectivity with disconnected boundary condi- 
tions 

We will now aim at transferring our results to the case of "disconnected boundaries," i.e. 
where no connections between the left-most and the right-most node are possible. This is 
the situation considered by Santi and Blough SBQ3], and one of our goals is to compare 
our results to theirs. We wish to show that the specific form of boundary conditions has 
no effect in the limit n —* oo, that is, our results from Sec. 13.2.31 hold for disconnected 
boundary conditions, too. 



3.3.1 Estimates 

Both models - the MANET with periodic and disconnected boundary conditions - are 
formulated on the same probability space, but are based on different events and random 
variables. For the disconnected boundaries, we consider the events k-DISCONN-DB, de- 
fined as 

Mi-DiscoNN-DB = {x E [0, 1]™ | x' i+1 — x • > p for exactly k values of i £ {1, . . . , n — 1}} 

(3.55) 

where x[ < x' 2 < ... are the sorted coordinates (a;*); as discussed, these events differ 
from k-DISCONN-PB only by the handling of the nodes on the boundary. The event 
k-DISCONN-DB, or rather its characteristic function, is certainly scaling and also sym- 
metric, since it only refers to the sorted coordinates. However, it is no longer translation 
invariant. Thus we can apply Eq. ()3.8|) . but no longer the results of Sec. 13.2.11 For exam- 
ple, we might calculate the probability of connectedness, i.e. of the event CONN-DB = 
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0-DISCONN-DB, as 

rl /■min(l,xi+p) /•min(l,a; n _i+p) 

-Pconn-db = n\ dxi dx 2 ... dx n . (3.56) 

Jo J xi J x n -i 

This integral can in principle be solved (for fixed n) and gives a piecewise-defined polyno- 
mial in p of degree n. (One might e.g. use computer algebra to solve it for realistic n.) 
However, a closed solution for arbitrary n seems to be out of reach. Instead, we will restrict 
to estimates of the difference between the periodic and disconnected boundary conditions. 

It is obvious that -Mk-DiscoNN-PB C -Mk-DiscoNN-DB j however, the opposite inclusion is 
not true: A configuration x G M k _ DISCO NN-DB might be disconnected at k + 1 places with 
respect to periodic boundaries, namely if the left-most point in x is not connected to the 
right-most one "via the boundary" . More precisely, let 

S := {x G [0, 1]™ | x[ + 1 — x' n > p, where x[ = min{xj}, x' n = maxjxj}}; (3.57) 

this is the event that the network is disconnected "at the boundary" . Then it is clear that 

^k-DISCONN-DB = A^k-DISCONN-PB U (^(k+l)-DISCONN-PB nS), (3.58) 

where the union is disjoint. This gives us the following inequality: 

-Pk-DISCONN-PB < -Pk-DISCONN-DB < -Pk-DISCONN-PB + -P(^(k+l)-DISCONN-PB H S) . (3.59) 

As an estimate, we can certainly say that 

-Pk-DISCONN-PB < -Pk-DISCONN-DB < -Pk-DISCONN-PB + P(S); (3.60) 

it then only remains to calculate P{S). 

The event S is most conveniently described in sorted coordinates; in fact, if x\ < . . . < 
x n , we can write its characteristic function as 

Xs(xi, ...,x n ) = Q{xi + (1 - p) - x n ). (3.61) 
Using Eq. ()3.8|) to calculate the expectation value, we get 

Ps = n\ / dx% I dx2 . . ■ / dx n 9{x\ + (1 — p) — x n ). (3.62) 

^0 J xi Jx n -i 

Let us split the integral in a sum Ps = I1 + I2, where I± covers the integration domain (p, 1) 
in the variable x±, and I2 covers the interval over (0,p). If x% G (p, 1), then the argument 
of the theta function is always positive, thus 

Ii — n! / dx\ I dx 2 . . . / dx n 1. (3.63) 

J p J xi J x n -\ 

Introducing new variables Z\ = (xi — p)/(l — p), and Zi — (1 — Xj) /(l — p) for i — 2, . . . , n, 
this reads 

I 1 = n\(l-p) n [ dz x I 1 dz 2 ...[ ' dz n = n\ (1 -p) n vol(K.)- (3.64) 

Jo Jo Jo 
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The volume of the n- dimensional standard simplex is known from Eq. (|A.15J) : our result 
thus is l\ = (1 — p) n . Now for the second integral, namely 



h = n\ / dx\ / dx2 ■ ■ ■ / dx n 9(x\ + (1 — p) — x n ). 

Jo J x\ J x n —i 



(3.65) 



Here we introduce new variables Z{ = (xi+x — x{) /(l — p) for i — 1, . . . , n — 1; this leads us 
to 

ro r 1-3:1 i- 1 ~ Xl z i n _l 

j 2 = n! (1 - p) n / dxi / d*i / dz 2 ... dz n -i 9(1 - } j z-). 



o 



(3.66) 

Note that the # function restricts the domain of integration for z to the (n — l)-dimensional 
standard simplex, which is completely covered by the integration since (1 — X\)/(l — p) > 1. 
Thus 

I 2 = n! (1 - p) 71 " 1 p vol(K-i) =np(l- p) n -\ (3.67) 
Combining the results for I\ and I 2 in Eq. ()3.6()j) . our result is: 



Lemma 3.10. Let k G N . T/ien 

Pk-DISCONN-PB < Pk-DISCONN-DB < Pk-DISCONN-PB + (1 — P)™ + n P (1 — P)™ ■ 

This estimate is certainly not very strict and might be improved, but it is already 
sufficient for our purposes: Note that in the limit n —>■ oo and p —>■ 0, we have ln((l — p) n ) = 
— np + 0(np 2 ), and ln(np(l — p) n_1 ) = ln(np) — np + O(p) + 0(np 2 ); thus we see from 
Eq. ()3.32|) that the difference between upper bounds and lower bounds in Lemma 13.101 
vanishes as np — Inn — > 77. This means that we can directly transfer the results from 
Theorem 13.81 to the case of disconnected boundary conditions. It is also straightforward 
to transfer the results for np — Inn — > ±00 from Theorem 13.91 Let us summarize this as a 
separate statement. 

Theorem 3.11. Let (p n ) be a sequence in M. + , and suppose there is an 77 e R suc/i £/iat 

np n — In n > n. 

Then, we have for every fc 6 N : 



p (n,pn) n^oo e / -m 

k-DISCONN-DB r,j e P 



In particular, 
In the case 
one has 



p {n,p n ) n^co (_ - V \ 

CONN-DB ^PV e P 

, n— ->oo , n— >oo 

np n — Inn ► +00 or np n — Inn > —00, 



p conn-db 1 or > respectively, Pcqnn-db °- 



Overall, this makes our claim precise that the choice of boundary conditions does not 
play a role in the limit n — > 00. 
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3.3.2 Comparison with the literature 

Now that we have established our results for the system with disconnected boundary 
conditions, we are in the position to compare them with existing results in the literature - in 
particular with those of Santi et al. [SBV01, SB02, SB03J who investigated the probability 
of connectedness using the same mathematical model, obtaining analytical estimates (with 
different techniques than ours) and also numerical results. 

We start with the analytical results. For comparison purposes, let us first state the 
following special case of Theorem 13.111 (We return to the parameters r and £ in place of 
p = r/£ here.) 

Corollary 3.12. Consider the 1- dimensional MANET model with parameters n, r = r(n), 
£ = £(n) and disconnected boundary conditions. If there is an e > such that for large n, 

nr > (1 + e) £lnn, then Pconn-db "~ > °° > 1. 

If, on the other hand, 

nr < (1 — e) £ Inn, then Pconn-db n ^°° > 0. 

For the same situation, Santi and Bio ugh [SB031 Theorem 7] state that, when expressed 
in our notation, 

(a) if nr = k £ln£ with some k > 2, then Pconn-db — > 1, 

(b) if nr = 2 £\ia£ and r(n) — > oo, then also Pconn-db — * 1 5 

(c) if nr = k £\n£ with k < (1— e) for some < e < 1 and r e 6(£ e ), then Pconn-db ~h 1> 

(d) if (nr)/(£\n£) -> 0, then P C onn-db t 4 1, 

in the limit n — > oo and £ — ■> oo, under the additional condition that r/£ — > 0. While 
these results are formally quite similar to ours, they are only compatible with them when 
the factor \n£ does not differ significantly from Inn. In practice, one will usually consider 
the case r = const., and here in fact no difference arises. However, from a mathematical 
standpoint, this is not implied by the conditions of the theorem; in fact, one may construct 
cases where the predictions of Corollary 13 . 1 21 are in conflict with the results from [$B03j. 
For example, consider the case n = k £ k+1 , r = £~ k In £, where k > 2. Then Pconn-db - * 1 
according to |(a)| but Pconn-db - ► according to Corollary 13.121 On the other hand, let 
£ = e n , r = e n /lnn. Then Pconn-db — > 1 according to Corollary 13.121 in contradiction 

to 

The present author claims that these differences are due to the fact that the argu- 
ments presented in jSBVOH ETRH21 [SB03J are inconclusive. The reader will find a detailed 
discussion hereof in Appendix iBl 

Another analytical result for connectedness was obtained by Piret |Pir91j in a similar 
situation: He modeled the nodes of a 1-dimensional MANET by a Poisson process of 
constant density d, and proved that for the radio range set to 



(3.68) 
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Figure 3.3: Numerical data for Pconn-db compared with the analytical results. Parameters 
are set to n = ^/i and r = ky£log 2 i, with different values of k. The numerical data was 
taken from jSBVOll Fig. 3]. 



with a constant k > 0, one has 

' 1 if k > 2 



l— >oo 

Connectedness * 



(3.69) 
ifife<2. 



Due to the use of a Poisson process, the node numner n is a random variable in this case 
and not a fixed number; however, for large i (and hence n) one would expect that n assumes 
its mean value id with low variance. In fact, setting n = id in Eq. ()3.68|) . the result (j3.69|) 
is just what Corollary 13.121 amounts to. 

Santi and Blough also presented numerical data for Pconn-db, derived by statistical 
simulation. In fact, their numerical method amounts to a Monte Carlo approximation of 
the integral 

-Pconn-db = / d n x Xconn-db(^)- (3.70) 

J[0,l] n 

We will use their results from [SBV01, Fig. 3] for comparison with our analytical results. 
(The data given in 1SB031 Fig. 2] is not much suited for our comparison, since most 
data points correspond to very low values of n or p.) Figure l3~3l shows these data series 
together with the upper and lower bounds from Lemma 13.101 where -Pconn-pb is given by 
Theorem 13.71 The numerical data is compatible with the analytical bounds within the 
level of precision that would be expected from a Monte Carlo type of approximation; since 
the underlying mathematical model is identical in both cases, any differences can only be 
due to numerical precision. 

In fact, Fig. 13.31 suggests that the exact values of -Pconn-db are much nearer to our 
upper bounds than to the lower bound Pconn-pb- This can be understood at least on a 
heuristic level: In a situation where the largest part of Q n (more than 80% of probability) 
corresponds to connected networks, it would be expected that most of the remaining parts 
of Q n fall into iWi_ DI gcc>NN-PB, and that M 2 _disconn-pb and higher disconnected events can 
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rather be neglected. (Note that according to Theorem IH.8[ the -Pk-DiscoNN-PB asymptoti- 
cally follow a Poisson distribution.) Thus, the estimate in Eq. (|3.60j) . where we replaced 
P(M 1 _ DISCO nn-pb H S) with P(S), is quite tight in this case. 
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Chapter 4 
Quality Measures 



When concerned with the structure of a MANET on a low level, i.e. related to the mere 
connectivity of nodes, the question arises by which quantitative properties the quality of 
the network should be described. Several measures have been proposed in the literature to 
that end. Most commonly, it is required that the MANET should be (strongly) connected 
with high probability; however, this requirement turns out to be quite strong, and so one 
may want to consider more general, in particular weaker, measures. 

In this chapter, we will introduce a general notion of so-called quality parameters for 
MANETs, and show that detailed results for specific parameters can be obtained at least 
in simple models. In particular, this allows us to discuss the scalability of such quality 
measures for large systems. 

We will first define more exactly what a quality parameter is in our context, and 
introduce a classification of such parameters according to their scaling behaviour; this is 
done in Sec. 14.11 We then consider several specific quality parameters in Sec. 14.21 and 
calculate their expectation value in the 1-dimensional MANET model. In Sec. 14.31 we 
will compare our results to numerical simulations conducted by Roth [Rot03j in a similar 
model. Lastly, in Sec. 14. 4[ give some examples of quantitative predictions for MANET 
design that follow from our analysis. 

4.1 General properties of quality parameters 

A quality parameter for MANETs in our context is a random variable Q : Q n — > M, or 
rather a family of such random variables (for different parameter values). The average 
"quality" of the MANET is then described by its expectation value Q = Ei[Q]. We will 
usually choose the range of Q to be [0, 1]; however, this is only a matter of convention. 

The definition of specific quality parameters naturally is very dependent on the usage 
scenario and application. However, there is one overall property that we wish to discuss 
in a general context: It relates to the scaling behaviour of the system, since we are usually 
interested in the limit of large MANETs (n — > oo). 

Let us consider the 1-dimensional MANET model from Chap. H3 for concreteness. If 
the quality parameter Q is an "intrinsic property" of the system, that is related to its 
behaviour in the bulk, then one might expect the following: If we take, say, two MANETs 
with identical parameters n, r, £, and couple them together - i.e., we join the two intervals 
and consider them as a single network with the double node number, allowing connections 
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between the two parts -, and if the original MANETs had a quality of Q = E[Q], then 
the joint MANET should have the same quality value Q, at least approximately for large 
systems. This would mean 

E[Q (n ' rA ] « E[Q {2n ' r ' 2£) ] or, equivalent^, E[Q M ] « E[Q {2n ' p,2 \ (4.1) 

referring to the normalized radio range. Of course, the same heuristic argumentation should 
hold when tripling the system size, dividing it into parts, etc.; more generally, the quality 
value should depend on np only, rather than on n and p independently. Let us formulate 
this more precisely. 

Definition 4.1. In the model of a 1- dimensional MANET, a family of random variables 
Q( n 'p) i s called intensive 1 if there exists a function Q : R + — > K with the following proper- 
ties: 

• Q is not globally constant; 

• Given v e IR + and a sequence (p n ) in R + such that np n — > v as n — > oo ; one has 

E[Q {n ' pn) ] -> Q{v). 

Here the first condition is introduced in order to exclude "trivial" intensive parameters, 
such as those where always Q —>■ when np — > const. Note that the parameter np = nrjl 
can be interpreted as the "non-statistical degree of coverage" of the network: E.g. np = 1 
means that the radio range of all nodes combined covers the interval [0,£] exactly once. 

By the above definition, we do not mean to say that only intensive quality parameters 
are relevant for our system, or that non-intensive parameters are not meaningful. In fact, 
such non-intensive quality parameters may be required for some applications. However, 
one should keep in mind that these parameters may not scale well for large systems: For 
example, if we need np — > oo in order to keep the quality level of the system constant as 
n — > oo, then this means that the average number of nodes per interval of length r needs to 
grow arbitrarily in the limit; thus we are likely to run out of local channel capacity. Hence 
applications which rely on a high quality level with respect to non-intensive parameters 
may not be feasible in networks with a high node number. 

4.2 Specific quality parameters 

We will now investigate a number of specific quality parameters and calculate their ex- 
pectation value in the 1-dimensional MANET model introduced in Chap. EJ where we 
will always refer to the case of periodic boundary conditions. Our choice of quality pa- 
rameters mainly follows a discussion by Roth [Rot03j , who introduced four such measures 
(segmentation, area coverage, vulnerability, and reachability) in the context of a numerical 
simulation. 

lr Thc usage of the word intensive is motivated by an analogy to statistical physics: Here a thermody- 
namic variable, e.g. a state parameter for a gas, is called intensive if it does not change when the system 
is divided into parts; examples include temperature, pressure, and particle density. 
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4.2.1 Connectedness 

One obvious choice for a quality parameter is the probability that the network is connected, 
which we had already investigated in Chap. El So, more formally, we set Q connectedness = 
Xconn-pb, where we know from Theorems 13.81 and 13.91 that 

E[Q C onnectedness] = ^CONN-PB ~> exp(-e^) as Tip - In 71 -> T/, 

^ as rip -> i/. (4.2) 

Thus Q connectedness is noi an intensive parameter. As discussed above, this means that 
applications relying on connectedness of the network will not scale well in large systems. 

Closely related to connectedness is the quality measure of coveredness, investigated by 
Piret |Pir91j in 1-dimensional systems. Coveredness (not to be confused with the area 
coverage parameter that we will discuss in the next section) measures whether each point 
in the interval [0,£] is covered by the range of at least one MANET node. It is clear that 
we need precisely yi < 2p for each next-neighbour distance yi to achieve that the interval 
is completely covered, while the criterion for connectedness is yi < p. Thus, coveredness is 
related to connectivity by 

(j( n <P) _ £)K 2 P) (A Q\ 

Coveredness ^Connectedness ' V / 

and we can apply the above result (|4.2|) accordingly. 

4.2.2 Area coverage 

Area coverage is the area A covered covered by the range of at least one MANET node, 
divided by the total area A to tai of the system: 

Q Coverage ~j • (^-4) 

A total 

Its expectation value may be understood as the probability that an external network node, 
with its position randomly chosen, will be able to connect to at least one of the n nodes 
of the MANET. 

In our 1-dimensional model, "area" is to be understood as the length of the corre- 
sponding line segments. Note that through dividing by A tota i = £, our parameter Qcoverage 
is scaling in the sense of Definition 13. 1| thus we may again pass to the normalized radio 
range and set A tota i =1. It is also easy to express Qcoverage in terms of next-neighbour 
variables: The distance yi leaves an area uncovered if yi > 2p; if so, the length of that area 
is yi — 2 p. Thus we get the following expression for Qcoverage : 

n 

^Coverage = 1 ~ ^{Vi ~ 2p) %i ~2p). (4.5) 

i=l 

In order to determine its expectation value, we will calculate 

E[{ Vi - 2p) 0{ yi - 2p)\ = [ dp T n '^{y) ( yi - 2p) 9(y t - 2p) (4.6) 

J [0,1]" 
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for each fixed i, where we will assume p < 1/2 (otherwise, we trivially have Qcoverage = 1)- 
Lemma IA.4I and Proposition IA.3I of Appendix IA.1I then yield 

E[( W - 2p) 6{ Vi - 2p)\ = (1 - 2p) n - 1 [ dp T n -°\y) (1 - 2p) Vi = -(1 - 2p) n (4.7) 

•/[0,1]" n 

Inserting into the expectation value of (|4.5|) . we obtain 

E[Q C overage] = 1 - (1 - 2p) n . (4.8) 

(Again, this is valid for p < |.) Using Taylor approximation ln(l — x) = — x + 0(x 2 ), we 
have 

In (1 - 2p) n = -2np + 0{np 2 ); (4.9) 

thus, in the limit np — > v (where n — > oo, p — » 0, and np 2 — > 0), the area coverage 
converges to 

E[Q Coverage ] -> 1 - e" 2 ". (4.10) 
This means that the area coverage is an intensive quality parameter. 



4.2.3 Segmentation 

The segmentation of a MANET counts the number of disconnected segments in the net- 
work, i.e. the number of subgraphs into which the network graph is separated: We set 

# of network segments 

Q Segmentation 7, 7 7 '-, ] • (4-1 1) 

of network nodes 

In order to take account of the periodic boundary conditions, we will count the strongly 
connected situation (the event CONN-DB) as having network segments. (This explains 
the slightly modified setting in Eq. 1)4.11)1 when compared with the original definition by 
Roth Rot03]], who defined 

# of network segments — 1 . . 

^Segmentation 7, 7 ' \ ] Z • \ 

# of network nodes — 1 

This difference is rather a matter of convenience and should not play a role in the limit of 
large systems.) 

Within our 1-dimensional system, it is easy to derive an explicit expression for the 
segmentation: We know that the event k-DISCONN-PB corresponds to a situation with 
exactly k network segments. Since these events are disjoint, and since their union (over 
k = . . . n) exhausts the sample space Q n , it follows that 

1 n 

^Segmentation = ~ ^Xk-DISCONN-PB (4-13) 

k=0 

and consequently 

1 n 

E[Q Segmentation] = _ ^^k-DISCONN-PB • (4-14) 
fc=0 
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The probabilities under the sum are known from Theorem IH.7I 

n [1/ P ] 



E[Q S egmentation] = " £ £(-1)'"** K ) T (1 "JP)^- (4-15) 

n k=o j=k \ / \3 / 

Now observe that in the sum over j, we may as well replace the lower limit with 0, since the 
binomial coefficient ({) vanishes for j < k. We may then exchange the order of summation 
and get 

E[Q S egmentation] = - £("1) J ( ) C 1 ~ 3 PT"' £("!)** [ i) " ( 4 ' 16 ) 

Likewise, we may replace n with j in the upper limit of the sum over k, since the summand 
vanishes for k > j as well as for j > n due to the binomial factors. Referring to Lemma [A. 91 
in Appendix IA.41 we know that 

E(-i)**(T) = { -1 ifj = 1, ( 4 - 17 ) 

^ \kj otherwise. 



k=0 



So in Eq. (|4.1fij) . only the summand for j = 1 remains. Assuming p < 1, that leads to the 
result 

E[Q Segmentation] = (1 ~ P)" ■ (4-18) 

With arguments as in Eq. ()4.9)1 . this means that in the limit — ► z/, 

EfQscgmcntation] ~^ e " 5 (4.19) 

so ^segmentation is an intensive parameter as well. 



4.2.4 Vulnerability 

The next quality parameter we will consider is related to the question how much the 
network quality or topology changes when a single node is removed from the network. 
Specifically, we define the importance of the network node with number j as 

Ij := max{0, (# segments with node j removed) — (# segments)}; (4.20) 

i.e. Ij is the number of network segments which are created by switching off node j in the 
current configuration. Nodes with Ij > make the network "vulnerable" against changes. 
This motivates to define the vulnerability of the network as 



Vulnerability 




(4.21) 



In our 1-dimensional model, the importance of a node is either 1 (if removing the nodes 
splits the respective network segment in two) or 0. The ordering of nodes is not of relevance 
for Eq. 1)4.21)1 : so we may describe the event j'-IMPORTANT (meaning that Ij = 1) directly 
in next-neighbour coordinates as 



^-.important = {V I (Vj-i < P) A (yj < p) A + yj > p)}, 



(4.22) 
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where the coordinate indices are understood "modulo n," i.e. y$ is identified with y n . We 
will assume n > 2 in the following, so that yj and yj~i are independent coordinates. Taking 
the complement of the set above, we can say that 

^'-IMPORTANT = 1 - ^(^-important) = 1 - P{Vj-\ > P V Vj > p V + yj < p). (4.23) 

On the last expression, we apply the inclusion-exclusion formula from Appendix IA.2j this 
yields 2 

^-IMPORTANT = 1 - P(.Vj-l >P)~ P(Vj > P) ~ P(Vj-l + Vj < P) + P(Vj >PA > p) 

+ P{Vj-i > P A j/j-i + % < p) + > P A + < p) 

- > P A > p A + < p). (4.24) 

The last three summands of this expression obviously vanish. Moreover, we know from 
Lemma 13.61 that 

P{Vj >P) = P(Vj-i >p) = (l~ PY~\ (4-25) 
PiVi-i > P A yj > p) = (1 - 2P)"" 1 ; (4.26) 

here we have assumed p < 1/2. Further, Lemma f A. 51 in Appendix IA. II shows that 



p (vj-i + Vi < p) = / d pl eq (y) °(p - Vj-i - Vj) 

J\0,l] n 



'[0,1]* 

= l-(l-pr 2 (l + (n-2)p). (4.27) 
Combining Eqs. ()4.24|) to ()4.27|) . we have shown that 

^-important = (np - 1)(1 - p) n - 2 + (1 - 2p)"-\ (4.28) 
Inserting into Eq. (|4.21jl . we have obtained that for n > 2 and p < 1/2: 

1 n 

E[Qvulnerabilit y ] = ~ ^"-IMPORTANT = (np - l)(l - p)^ 2 + (1 - 2p) n ~ 1 . (4.29) 
U 3=1 

A Taylor approximation (as in the previous sections) then leads us to the following asymp- 
totic behaviour in the limit np — > v: 

E[Q V ulncrability] ~> (f ~ l)e~ V + C~ 2v . (4.30) 

Thus, the vulnerability is an intensive quality parameter as well. 
4.2.5 Reachability 

The reachability parameter is concerned with the number of nodes that can be reached 
from a given node (in a multi-hop fashion), or, alternatively speaking, with the size of the 
segments of the network. We define the reachability of some fixed node j as 

„ # of nodes reachable from node j 

Ra := — -. 4.31 

n 



2 More specifically, we apply Theorem lA.6l with respect to the event C\ and for n — 3 (with notation as 
in the theorem). 
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Here we do not count the node itself as reachable, unless the network is strongly connected 
(i.e. the node can "reach itself" via the boundary). We define our quality parameter, the 
average reachability, as 

1 n 

^Reachability = _ ^^.Rj- (4.32) 

U 3=1 

Again, we have introduced a slight difference compared to the original definition by Roth 
[Rot03j which accounts for the periodic boundary conditions and vanishes for n — > oo. 
Following our above discussion, the value of (^Reachability is 

• 1 in the event CONN-PB, 

• (n - l)/n in the event 1-DISCONN-PB, 

• more generally, n~ 2 52*=i ~~ -0 in ^ ne even t k-DISCONN-PB, k > 1, where 6, are 
the sizes of the k network segments. 

To get a more explicit description of the latter case for k > 2, we define the events 
SEGMENT-m-b, where m G {1, . . . , n}, 6 G {1, . . . , n — 1}, which describe that a segment 
of the network begins exactly at node m, extending "to the right," and has a size of exactly 
b nodes. (The node indices are counted in sorted coordinates, and are defined modulo n.) 
This can be formally expressed as 

m+b-2 

XSEGMENT-m-b (?/) = %m-l ~ p) %m+b-l ~ p) \ [ 9(p ~ Vi) ■ (4.33) 

i=m 

It is then easy to sum over the size of the segments: Since the events SEGMENT-m-b are 
obviously disjoint from CONN-PB and 1-DISCONN-PB, one simply has 

n — 1 ^.^--kb(b—l) (a?a\ 

^Reachability — XCONN-PB H Xl-DISCONN-PB + 2 XSEGMENT-m-b • (4.34) 

m=l 6=1 

Since the expectation value of the first two summands has already been calculated in 
Chap. El it only remains to calculate -PsEGMENT-m-b in order to determine E[Q Reachability ]. 
Using the definition in Eq. ([4.33)1 . and applying Lemma lA.41 twice, we see that for n > 2 
and p < 1/2, 

m+b-2 



[ d^{ y )e{v m ^-p)e{ Vm+b ^-p) TT e{ P - l 

J m n i=m 

„ m+b—2 
„ m+b-2 



[l-pT-\l 



=(l-2p)»-! 
\n-l 



(1 - 2p) n - 1 P(y m < p A . . . A y m+b ^ 2 < p'), (4.35) 
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where p' — p/(l — 2p). For determining the probabilities P(y m < p' A . . .), we once again 
use the inclusion-exclusion formula 3 of Appendix IA.2I 



P(y m < p A ... A y m +b~2 < p) = 1 - P(y m > p' V . . . V y m +b-2 > p) 

6-1 / . -i \ 6-1 



1=1 V / ,_n 



j=l \ / j= 

where 

Sj= Yl P(y mi >p' A...Ay m . >p'). (4.37) 

{mi ,...,m,j }c{m,...,m+6— 2} 

We already know the probability under the sum by Lemma T3. 61 Applying this result leads 
us to 

[i/p'l 

"' /] jp'T- 1 - (4-38) 

3=0 

Now we can assemble our results, together with the expressions for Pconn-pb and -Pi-disconn-pb 
from Theorem 13. 7^ in order to determine the expectation value of Eq. (j4.34jl : This gives 

n — 1 n n l bib — 1) 

E[QReachability] =PcONN-PB H Pl-DISCONN-PB + ^ PsEGMENT-m-b 

m=l 6=1 

[!//>] ,_ N _ 1 [VP] 

\n— 1 



P(ym<p'A...Ay m+b -2<p')= 52(-lY[ ■ )(1 



Et-^Qti-^r'+^Et-iwQti-^)' 

+ - 2,)-' £ ( 6 " M (1 - (4.39) 



6=1 j=0 

where p' — p/(l — 2p), and we assume n > 2, p < 1/2. 

While this explicit expression is rather complicated, we can derive a much simpler result 
for the limit n — > oo, where we consider np — Inn — > 77 as in Sec. 14.2.11 We already know 
the limit values of Pconn-pb and -Pi-disconn-pb from Theorem 13.81 It is also easy to see 
that 

ln(n 2 (l - 2P)"- 1 ) = 2 Inn - 2np + 0{p) + 0{np 2 ) -»■ -2r/, (4.40) 

so the factor n 2 (l — 2p) n ~ l converges to e~ 2r? . It remains to determine the asymptotic 
behaviour of the sum over b. The idea here is to understand the sum (for large n) as the 
approximation of a Riemann integral, where the integration variable j3 = b/n ranges from 
to 1. Since the calculation is somewhat involved, we state it as a separate lemma. 

Lemma 4.2. Let rj G R, (p n ) C M + such that np n — Inn — > 77 as n — > 00, and Ze£ 
Pn := Pn/(1 — 2p„). Tnen one has 



6=1 j=0 \ J / JO 



3 More precisely, we use Theorem I A . 61 with respect to the event C\ and with (6 — 1) in the place of • 



4.2 Specific quality parameters 



41 



Proof. In the following, we keep rj fixed and set 

/(/3)=exp(-/3e-"), (4.41) 

up) = E " ^(i -i^r 1 , (4.42) 

i=o 

and o B =j!^ 1 Vl- J y n ) B - 1 . (4.43) 

We obviously have |/(/3)| < 1 for (3 G [0, 1], and we also know that \ f n (P)\ < 1 for /? = b/n, 
b G {1, . . . , n}, since the f n {b/n) are defined as probabilities (cf. Eq. ()4.38|) : we can easily 
extend this to the case b = n). This is useful for simplifying the proposition of the lemma: 
Since 

'2^/.(UI<--0 (4.44) 



rr n 



and 

n i2 



E 6 ^'W ll^E-S-^O, (4.45) 

6=1 

we can equivalently prove that 



rr rr * — ' n n 

6=1 6=1 



' n J v rr / n 

6=1 1/0 



- / d/3/3 2 f(l3)\ ^^0. (4.46) 
However, since / is integrable, it is clear by the definition of the Riemann integral that 

IElsf/0-/ ^/? 2 /(«|^0. (4.47) 

6=1 J0 

Thus, it only remains to verify that 

£%«(-)-/(-) 1^0. (4.48) 
^— ' rr 1 n' n 1 

6=1 

To that end, we need an estimate of \f n (b/n) — f(b/n)\ that is uniform in b. We will 
construct this estimate by refining the methods developed in the proof of Theorem 13 .8| 
using notation as introduced there. 4 

Regarding the terms aj, we can certainly say that for j < b — 1, 

'( > j') = Mli - li "' F -''' (4 ' 49) 

independent of 6; the same is true for j > 6 — 1 (where the binomial coefficient vanishes). 
We can then apply the same construction that lead to Eq. ()3.43|) . Thus, for given e > 0, 
we can find j and n such that for any n > n , 

[l/pn] 

\J2—r a ^ 2e - ( 4 - 5 °) 

3=30 



J Note that the parameter k in Theorem 13.81 must be set to for our purposes. 
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and at the same time, for any 6 E {1, . . . , n}, 

i-iy.b 



E^(r")'|S, (4.51) 



J=JO 



(Note that we can find such an estimate independent of 6, since the power series Y2j x ^ I ft 
converges uniformly on the interval [— e~ v ,0}.) 

Now it remains to handle the terms for j < jo; we have to find a uniform estimate for 
\(-e~ r ') : ' — aj \ for all b at fixed j. Let us first consider those terms where b > en, where we 
can assume that j < en (possibly after increasing n ). We know that 

a 3 /( b -e~y = n'e*(l - jX)- 1 . (4.52) 

Only the first factors in this expression depend on 6; they are 

(6-1)! 1 6-1 b-j 



(b-l-j)\b> 6 "' 6 
Each of the factors of the form (6 — i)/b converges to 1, more explicitly: 



(4.53) 



1^-11 = 1^- («4) 
6 6 en 

Thus we can control the convergence of these factors independent of 6 (with jo still being 
fixed). Moreover, we find - just as in Eq. 1)3.49)) - that 

nV^l-jp;)- 1 ^!, (4.55) 

where the term does not depend on 6. Thus the convergence of aj/(^e~ v )i — > 1 is uniform 
in 6, given that 6 > en. Summarizing this with Eqs. ()4.50|) and ([4.51)1 . we have found that 

Ve > 3m \/n > m V6 G {[en] + l,...,n} : |/ n (-) - /(-)| < 4e. (4.56) 



n n 



For 6 < en, we will use the rough estimate 



|/»(")-/(-)|<2. (4.57) 
Now combining these bounds, we can establish Eq. ()4.48|) : For n > ni, we have 

n 1.2 i. 7 n 1.2 [ £n l 7,2 i i 

E^I/4)-/0l^e E ^^ E ^<44 n + 2l W <6, (4.58) 

6=1 b=[en]+l 6=1 

This finally proves Eq. ()4.48|) and hence the lemma. □ 

Of course, the integral that we established as a limit value in the above lemma is easy 
to solve (twice integrating by parts): One has 

i 1 
dp P 2 exrj(-Pe~ v ) = - exp(-/?e- ,^ )(2e 3, ' + 2e 2 "/3 + e^? 2 ) 

= - expt-e^Xe" + 2e 2r} + 2e 3?? ) + 2e 3 ". (4.59) 
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parameter 


expect at 
at finite n > 2, p < 1/2 


ion value 

asymptotic 


intensive? 


Q Connectedness 


[i/p] / \ 

3=0 KJ J 


exp(— e _T? ) 

as np — Inn — > n 


no 


Q Coverage 


1 - (1 -2p) n 


1 - e~ 2v 
as np — > z/ 


yes 


Q Segmentation 


(1 - pf- 1 


as np — > z/ 


yes 


Q Vulnerability 


[np- 1)(1 - p) n ~ 2 
+ (1 — 2p) 


[v - l)e~ v + e- 2 ^ 
as np — > z/ 


yes 


^Reachability 


see Eq. (|OS|) 


2e^ - (1 + 2e r ')exp(-e- ? ') 
as np — Inn — > 77 


no 



Table 4.1: Overview of the results for quality parameters 



Now collecting our results on E[QR eac h a biiit y ] in Eq. (|4.39|) . where the limits for Pconn-pb 
and -Pi-disconn-pb are known from Theorem IH.8| we can establish that 



E[<5Rcachabiiit y ] -> ^ - (1 + 2e v ) exp(-e n ) as np - Inn -> r^. 



(4.60) 



By a monotony argument similar to the one which lead to Theorem 13. 9[ we can show that 
Ef^Reachabiiity] — ► as np — >• z/; so the reachability is noi intensive. 



4.3 Comparison with simulations 

We will now aim at comparing our results on quality parameters, which are summarized 
in Table l4~Tl to the simulation data obtained by Roth [Rot03j. 

In contrast to the quite simplistic assumptions of our model, Roth aimed at a more 
realistic network topology; he chose part of the map of the Downtown Minneapolis shopping 
center as the basis for his simulation (cf. Fig. 14. lj) . This shopping center consists of a 
number of towers which are connected on the first floor via bridges, so-called "Skyways" ; 
we consider users with wireless devices moving along these paths (see Fig. 14. 2|) . 

This model is in a way quite similar to ours and largely makes the same overall assump- 
tions: Network nodes move independently at random on 1-dimensional paths; the radio 
range of all nodes is equal with a sharp cutoff at radius r. However, there are a number of 
important differences: 

First, while we based our analysis on a static model (assuming ergodicity for mobile 
nodes), Roth considered an explicit motion model: Users move at constant speed along a 
line segment, and choose a new speed and direction once they have reached the end of a 
segment. Certainly, one would expect that this model also leads to an equal distribution 
of nodes on the line segments in the long run; however, this is not explicitly modelled. 

Second, Roth considered a 2-dimensional radio propagation, in contrast to our 1- 
dimensional model; i.e. two nodes are connected when their distance is smaller the r 
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Figure 4.2: Idealized map of the Skyways. Dashed line segments were not considered for 
determining the effective length I (see text). 



4.3 Comparison with simulations 



45 



on the plane rather than along the line segments. (No shielding by buildings, walls, etc. 
between the different paths was taken into account.) In most cases, this is equivalent to 
our 1-dimensional propagation, since neighbouring line segments are usually further than 
r apart (cf. Fig. 14. 2|) : however, there are some exceptions. We will discuss this in more 
detail below. 

Third, as already noted, the topology of the line segments is much more complex than 
in our simplistic model, including both open and closed curves. 

Before we can compare our results to those of Roth, we must first determine the param- 
eters of our model that correspond to the situation considered by Roth. The radio range 
was chosen as r = 30m (the indoor communication range of IEEE 802.11b Wireless LAN), 
which we can directly transfer to our situation. The system length I is more difficult to 
determine: While it might seem obvious to set I as the total length of all line segments in 
the system (see Fig. I4.2J1 . there are two corrections we wish to make. These are due to the 

2- dimensional propagation model used by Roth. 

On the one hand, Roth's model allows communication between nodes on parallel (or 
nearly parallel) line segments whose distance is less than the radio range. In our model, 
however, nodes can only communicate in direction of the line segment. Thus the range 
of a network node covers additional segment length in Roth's calculations, the more the 
nearer such parallel line segments are located. We will roughly accommodate this effect by 
the following procedure: Whenever two parallel line segments in the map are not further 
than r/2 apart, we will only count one of them for determining the total system length I. 
The line segments that were left out due to this procedure are marked as dashed lines in 
Fig.EOl 

On the other hand, there is another effect at those points were at least 3 line segments 
meet. Due to the 2-dimensional propagation model, nodes which are located near such 
a point can reach other nodes in line segments of approximately 3r in length (lr in each 
direction); in our model from Chap. however, nodes can only reach an "area" of 2r in 
length. In order to compensate this difference, we will subtract lr from the parameter i 
for each such point on the map. There are 30 points of the mentioned type on the map, 
not counting line segments that were left out due to the procedure described earlier. This 
leaves us with an effective length of 

i = 3363 m - 30 • 30 m = 2463 m. (4.61) 

Of course, these "ad hoc corrections" are only very rough and cannot be traced back 
directly to the statistical description. They also do not account for all effects that relate to 
differences between the models - for example, the 2-dimensional radio propagation certainly 
has an effect that relates to points where only 2 segments meet, while the effect around the 

3- segment points may have been over-estimated; also, we do not account for the increased 
density of nodes in the areas where two line segments run in parallel. However, we shall 
see that with the corrections introduced, we can already get a good match between the 
results that the two models predict. 

After having fixed the parameters, let us now turn to a direct comparison of the data. 
Roth did not consider connectedness as a quality parameter, since in fact (as discussed 
above) strong connectivity would be a quite strict condition for networks of reasonable 
size. So we will discuss coverage, segmentation, vulnerability, and reachability. For all 
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• simulation results 
analytical results 



Coverage 




100 200 300 400 500 600 



(a) linear 




100 200 300 400 500 600 



(b) logarithmic 



Figure 4.3: Analytical and simulation results for the area coverage parameter 
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• simulation results 
analytical results 



Segmentation 




100 200 300 400 500 600 



(b) logarithmic 



Figure 4.4: Analytical and simulation results for the segmentation parameter 
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analytical results 




n 



100 200 300 400 500 600 



(a) linear 

( ^Vulnerability ) 




n 



100 200 300 400 500 600 



(b) logarithmic 

Figure 4.5: Analytical and simulation results for the vulnerability parameter 
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Q 



Reachability 
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. 6 
. 4 
. 2 



• simulation results 
analytical results 



100 200 300 400 500 600 



Figure 4.6: Analytical and simulation results for the reachability parameter 



these parameters, we will compare the numerical results of |Rot03j with our explicit re- 
sults listed in Table 14.11 where we will use the exact formulas rather than the asymptotic 
approximations. (In most cases, the difference between the asymptotic approximation and 
exact value is however so small that it would hardly be visible in the graphs.) 

Let us start with the area coverage parameter, shown in Fig. 14 .HI The linear plot 
shows that both models nearly agree in absolute values for n = 50 and n = 100, and 
in the asymptotic behaviour as n — >• oo (where both graphs approach 1), while there is 
some difference at medium values of n. However, the logarithmic plot reveals that our 1- 
dimensional model systematically differs from Roth's simulation, which shows a much lower 
area coverage at high n. An explanation for this difference might be boundary effects in 
Roth's model: Possibly, some peripheral parts of the Skyways were not as densely covered 
with nodes as one would expect from the equal distribution. Still, the absolute difference 
between the models is below 5%, and the models agree with respect to their qualitative 
behaviour. 

The data for segmentation is shown in Fig. 14.41 It shows a good fit between the models, 
both on the linear and logarithmic scale. In particular, Q segmentation decays exponentially 
with n quite precisely, which is visible in the logarithmic plot; this is exactly the behaviour 
predicted by our simpler model. 

Figure 14.51 compares the data for Q vulnerability For this parameter, we also obtain a 
good fit between the two models across the range considered for n, except perhaps for the 
case of very few nodes (n = 50). 

The last parameter - reachability - is shown in Fig. 14.61 While the qualitative be- 
haviour agrees between the models also in this case, there are noticeable differences in 
the absolute value of ^Reachability: In the range of medium n, it seems that in the simple 
1-dimensional model, approximately 50-100 nodes more are needed to achieve the same 
reachability as in the simulation by Roth. This leads to absolute differences of up to 0.3 in 
^Reachability between the models. Taking into account that the average number of network 
segments agrees between the models (cf. Fig. l4.4J) . this points to the fact that at least some 
particularly large segments occurred in Roth's simulation that are not predicted by our 
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1-dimensional model. This is possibly explained by the fact that Roth's model allows com- 
munication between parallel paths; while we compensated this partially by counting only 
one contribution to I from two parallel paths, this still amounts to an increased density of 
nodes in those areas that would not correctly be described by an equal distribution. 

Certainly, it would be possible to gain a better and more quantitative understanding 
of the difference between the two models by repeating and modifying the simulations of 
[Rot02j, and by refining the construction in Chap. El and 0] in order to include more complex 
situations. However, such an analysis lies beyond the scope of the current work. 

In conclusion, it seems that the numerical results in |Rot03j can be reproduced in our 
more simple model at least in a qualitative sense, and in large parts also quantitatively. 
It should be emphasized that this does not amount to a comparison with experiment; we 
merely compared our results to a different mathematical model, which is partially based on 
the same simplifying assumptions (e.g. a homogeneous radio range for all nodes). Still, the 
material of this section may support the claim that the predictions of our 1-dimensional 
system are stable with respect to some changes in the modelling decisions. Differences with 
respect to details of the propagation model could be compensated by a simple change in 
the system parameters. 



4.4 Quantitative predictions 

More explicitely than the results known in the literature, our asymptotic approximations 
allow us to make quantitative predictions for the quality of 1-dimensional MANETs under 
the given modelling assumptions, or, more importantly, to find appropriate system param- 
eters required to reach a certain quality level. This section gives some examples to that 
end. 

Assume in the following that the length i of the MANET and the radio range r are 
given. We want to find the minimum node number n needed to obtain different quality 
levels, where we restrict our attention to the case of large MANETs; i.e. we will use the 
asymptotic formulas for quality parameters from Table 14.11 on page E31 

Let us start with connectedness. Given some required quality level Q connectedness, we can 
directly obtain the associated value rj by r\ = — In (— In Q connectedness)- It remains to find 
n such that rj = nr /£ — Inn. Given r/£, this solution needs to be calculated numerically, 
which is however easy to do (e.g. using Newton's algorithm). 

For area coverage and segmentation, the required value of v and hence of n = u/r is 
directly obtained from Qcoverage an d Q segmentation without further complications. For the 
vulnerability, we need a numerical inversion of z/ i — > (i/ — l)e~ u + e~ 2u in order to obtain n 
from ^Vulnerability (One usually obtains two such solutions for n - cf. Fig. 14.51 -. where we 
are interested in the greater one.) Likewise, for the reachability parameter, a numerical 
inversion of x i— > 2x — ( 1 + 2x)e~ 1 ^ x gives us the required value of x = e v ; we then proceed 
as above in order to calculate n from rj. 

All these calculations can be performed with standard techniques (Newton's method, 
regula falsi) and without excessive need for computing capacity. In fact, the evaluation 
would be feasible even on a mobile device with very limited CPU power, should this become 
necessary e.g. within a distributed algorithm. 

Table [Ol shows some numerical examples for a MANET of I = 1 km in length, using two 
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criterion 


minimal node number n 




IEEE 802.11 WLAN 


Bluetooth 




(r = 30 m) 


(r = 10 m) 


^% Connectedness U.Jy 


261 


906 


Qcovcra^e > 0.9 


39 


116 




77 


231 


^Vulnerability ^ 0.1 


102 


304 


Q Reachability > 0.9 


173 


650 


^Connectedness > 0.99 
v ^ v^onncciecuicss 


349 


1167 


Qcoveraee > 0.99 
coverage ^ • ^ ^ 


77 


231 


^!/"*crm o vi ^ i t i . 1 


154 


461 


Q Vulnerability < 0.01 


209 


627 


^Reachability > 0.99 


226 


804 



Table 4.2: Quantitative predictions for a 1-dimensional MANET (£ = 1000 m). 



different radio ranges (for IEEE 802.11 WLAN and Bluetooth radios) and various quality 
criteria. As expected, the non-intesive parameters (connectedness and reachability) lead 
to criteria that are particularly demanding in terms of node density. For example, if one 
requires 99% probability of connectedness in a Bluetooth-based MANET, then more than 
1.100 network nodes are needed, which is more than one node per meter of network length 
- a threshold that would probably be hard to reach in practice. 
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Further Directions 



This chapter discusses extensions of our results to more complex situations. To that end, 
Sec. 15.11 presents a variation of our 1-dimensional MANET model in which the network 
nodes may be switched off at random. Sec. l5.2l then gives a summary of the results obtained 
in the current work, as well as an outlook to higher-dimensional systems and the description 
of time dependence. 



5.1 A network with varying node number 

As a simple example of how our method can be generalized to more complex behaviour, 
let us consider the following situation: In the 1-dimensional MANET, we introduce a 
varying node number by allowing each network node to be switched off at random. This 
corresponds to a user turning off their device e.g. for power saving reasons. We will assume 
that at any fixed time, each device is switched on with probability p (where the devices 
are independent of each other). This is reflected in the model by adding a sample space 
^internal = {0, 1} for each node, where the value corresponds to the device being switched 
off. We thus consider the sample space 

0>,vn= ([0,1] x {0,1})". (5.1) 

We extend the probability measure by adding a discrete distribution for each of the addi- 
tional coordinates Z{ G fi internal (i = 1, • • • ,n)] the expectation value of a random variable 
Fvn : ^ n ,VN — * K then is 

In „ 

E[^vn]= J2 (T\P H ( 1 -P) 1 ~' i ) d n xF VN (x 1 ,...,x n ,z 1 ,...,z n ). (5.2) 



z 1 ,...,z„=0 i=l 



Following our motivation, we can define our random variables of interest (i.e. the quality 
parameters) quite easily: We want that for our quality measures, only those nodes with 
Zi — 1 are counted. That is, for a given family of random variables : O n -> 1 on the 
original MANET (with fixed node number), we define a variable Fvn on the new sample 
space fi n ,vN by 

F<J${x,z)=F< n '\y), (5.3) 
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where n' = J2i z i> an d V = (vi, ■ ■ ■ , Vn 1 ) lists those variables Xi for which Zi = 1. This 
definition is unambiguous if the are symmetric, which was the case for all our quality 
parameters. 

For this specific choice of random variable i*vNj the expectation value from Eq. (|5.2j) is 
somewhat simplified: We can integrate over all variables that do not appear in F^ n '\ and 
make use of the fact that F^ n > does not depend on the Zj. This leads us to 

E[F VN ] = V] p^ Zi (l-p) n ~^ Zi [ d n 'yF^\y). (5.4) 

Since only the sum of the Z{ is relevant in this expression, we can replace the multiple sum 
by a single sum over n'\ 

E[Fvn] = (^/(l-Pl^E^'']. (5.5) 

n'=0 ^ ' 

Clearly, one would expect that for large n, the MANET with varying node number will 
behave like the MANET with fixed node number, but at the parameter value pn in place 
of n. Mathematically, this is a consequence of the central limit theorem. We shall show 
this precisely at least for some parameters of interest. 

Theorem 5.1. Let Q^ n,f ^> be a family of random variables for the 1-dimensional MANET; 
assume that Q^ n ' p ^ is scaling, symmetric, and intensive with limit function Q. Moreover, 
let Q^'ti be bounded in the sense that there exists a constant M > such that 

\/n € N Vp e M + Vcj G tt n : \Q {n ' p \u)\ < M, 

and suppose that the convergence Q —>■ Q is uniform in the following sense: 

Ve > 3n 35 > Vn > n : \np - v\ < 5 => \E[Q {n ' p) ] - Q(u)\ < e. 

Let be the corresponding random variable for the MANET with varying node number. 

Then, for each sequence (p n ) with np n — ► v > 0, one has 

E[QvX°] -»• Qipv). 

Note: It can easily be shown that the above condition of uniformity is fulfilled in all 
our examples of intensive quality parameters. Also, all our parameters were bounded by 
definition. An analogous theorem for our non-intensive parameters might be stated, but 
we will not discuss this in detail. 

Proof. In view of Eq. (|5.5jl . our task is to show that 

E ( n V'(l -p) n - n '\n_Q {n '' Pn) \ - Q{pu)\ o. (5.6) 

n'=0 \ U ' 

To that end, let e > be given. Further, let A > (its value will be specified later). Let 
o~ n = ^Jnp{l — p), a n = [np — Xa n ) , and j3 n = [np + Xa n ] . Applying the de-Moivre-Laplace 
theorem (cf. Theorem IA.8l in Appendix IA.3|) . we know that 

lim E ( n )\p n '^ - p) n ~ n ' = $ (- A )- ( 5 - 7 ) 

n'=0 V 7 
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(See Eq. (IA.25J) for the definition of $.) Since with Q, also its limit function Q must be 
bounded, we can thus obtain for large n: 



n'=0 

Likewise, we see that 



[™\v n '^-v) n ~ n '\E[Q {n ' ,pn) ] ~Q(pv)\ < 2M($(-A) +e/M). 



(5i 



J2 (1 -p) n ' n '\^{Q {n '' pn) ] -Q{pu)\ < 2M($(-A) +e/M). (5.9) 

It remains to control the sum over n! e + 1, . . . , j3 n }. For these values of n', we certainly 
have 

\n'p n -pv\ < p n \n' - np\ + p\np n -v\< Xp n a n + p\np n - v\. (5.10) 

Since p n = 0(l/n), a n = 6(y / n), and np n — >■ i/, we can achieve that |n'p„ — pv\ < 5 
for sufficiently large n, where S is the value used in the uniformity assumption. This 
assumption then guarantees that |E[Q^ n ' ,Pn ^] — Q{pv)\ < e for large n and a n < n' < (3 n ; 
thus 

r f )p n '(l-pT- n '\nQ {n '' Pn) }-Q(pis)\ <e. (5.11) 

n'=a n +l ^ ' 

Combining Eqs. ()5.8|) . (|5.9|) . and f|5. 1 1|) . and choosing A large enough such that $(— A) < 
e/M, we have achieved the desired result. □ 

The above theorem says that for the model with varying node number, we can apply 
the results from Chap. 0] directly if we set the node number in those results to np (i.e. to 
its statistical mean). While this is not very surprising, it means that our model is stable 
(to some extent) against changes in the assumptions; we can accommodate the extra effect 
by merely modifying one of the system's parameters. 

For illustration, let us discuss the above findings in one concrete example, namely 
the segmentation parameter (^segmentation introduced in Sec. 14.2.31 Here we know from 
Eq. gHED that 

E[Q Segmentation] = (1 ~ p)" ■ (5.12) 

Inserting into Eq. ()5.5|) . we can explicitly calculate the segmentation for varying node 
number: 



E[Q S egmentation,VN] = £ ( "V^l " p)^' (1 ~ pf'' 
n'=0 ^ ' 



1 .(p(i- p ) + (i- P ))» = £-Eer (5.i3) 



l_ p V- r, V r„ 1 _ f) 

In the limit np — > u, it follows that 

E[Qsegmentation,VN] C PV = (^Segmentation {pv) , (5-14) 

as expected. 
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5.2 Conclusions and outlook 

In the course of the present work, we have analysed a 1-dimensional MANET system 
with statistical methods. Using a number of symmetries of the system, the mathematical 
description of connectivity properties could be much simplified. It turned out that the 
model was explicitly solvable when boundary effects were neglected (through the use of 
periodic boundary conditions). In particular, we were able to obtain an explicit expression 
for the probability of connectedness for given parameters, and analyse this expression in 
the limit of large MANET size. This improves the results known in the literature for 
1-dimensional systems. 

We then analysed a number of different quality measures for MANETs. In general, 
quality parameters could be classified into intensive parameters (with good scaling prop- 
erties) and non-intensive ones (which possibly lead to scalability problems). We were able 
to obtain explicit results for all of the parameters in the simple 1-dimensional model. Our 
results agree with the numerical data known in the literature. 

Our results can serve both as a qualitative and quantitative guideline for the design 
of 1-dimensional MANET systems, in particular for sensor networks. Due to our explicit 
results for the expectation value of quality parameters, it is easy to choose the radio range 
or node density in a MANET such that it reaches the desired quality level. In particular, 
this applies to the asymptotic formulas; they are certainly simple enough to even allow 
computation on the mobile devices themselves. 

Further, the methods we have developed should be applicable also to other quality 
parameters, in case they are desired for specific applications: As long as these parameters 
can reasonably be expressed in terms of the next-neighbour coordinates, it should be 
possible to apply the techniques of Chap. 0]in order to obtain their expectation value. 

Certainly, we have merely treated a small part of the problems and obstacles that 
may limit the quality and scalability of MANETs. In particular, we have not dealt with 
questions of routing, throughput, and all aspects explicitly related to mobility. Thus, our 
results should be regarded as a upper bound to MANET quality, in the sense that additional 
problems might be faced on higher layers. 

Our specific 1-dimensional model is quite simplistic in its assumptions, and it would 
certainly be worthwhile to study some extensions in order to explore the stability of our 
results against changes in the model. Apart from an inhomogeneous spatial distribution 
of the nodes, it would be particularly interesting to analyse nodes with a varying radio 
range, which might be caused e.g. by local interference, changes in antenna positions, or 
shielding. In analogy to Sec. 15. 1[ this could be modelled by introducing additional random 
parameters into the formalism which control e.g. the radio range between each pair of 
nodes, or only between next neighbours. Still, one would expect that under reasonable 
assumptions, the extended model could effectively been reduced to the known situation by 
application of the central limit theorem. 

It would also be desirable to extend our findings to 2-dimensional and, with some 
limitations, to 3-dimensional MANET systems. In fact, some of the results can easily 
be generalized: Let us consider the area coverage parameter Qcoverage- Assume that n 
nodes with circular radio range r are distributed equally (and independently) to the cube 
[0,£] d (d e N), considered with periodic boundary conditions. We can certainly say that 
E[l — Qcoverage] is the probability that a dedicated point, distributed at random to [0,£] d , 
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Figure 5.1: Long-range dependence of the vulnerability in 2-dimensional networks. The 
importance of the node marked with a solid arrow depends on the position of the node 
marked with a dotted arrow, and vice versa. 

will fall into the range of none of the MANET nodes. The probability for the dedicated 
point to fall into the range of one specific node, however, is simply Qp d , where p = r/£ as 
usual, and q is the volume of the unit sphere in d dimensions. Due to the independent 
distribution of network nodes, we obtain 

EfQcoverage] = 1 - (1 - C d p d ) n -> 1 - aS n P d -> U, (5.15) 

in generalization of our 1-dimensional result in Eq. ([4.10)1 : we have 

4 

Ci = 2, C 2 = 71, C 3 = -7T. (5.16) 

The results for other quality measures, in particular for connectedness, do not transfer 
that obviously however: Since the next-neighbour coordinates cannot be used in the same 
way in higher dimensions, we would first have to find appropriate new coordinates in order 
to transfer our methods. On the other hand, similar results would be expected to hold; 
cf. the numerical results by Santi and Blough [SB03 and the analytical estimates by 
Bettstetter |Bet02j . 

It should also be noted that certain quality measures somewhat change their nature in 
d > 1 dimensions: As an example, consider the vulnerability parameter (cf. Sec. 14.2.4)1 . 
In the 1-dimensional situation, the question whether a node is "important" for network 
connectivity is determined by its two associated next- neighbour distances, and hence we 
may say that it is a local property. In d > 2, however, it may happen that the importance 
of a node depends on the structure of the network at a very remote place (see Fig. 15.1)1 . 
Since we want our quality measures to reflect the behaviour in the bulk network, it might 
even be necessary to change the definition of the quality parameters in higher dimensions. 

Up to now, we have only considered static deployments of network nodes, taking account 
for mobility only through our assumption of ergodicity. While for the quality parameters 
we considered, we are in good agreement with simulations that rely on an explicit motion 
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model (cf. Sec. 14.3)1 . our framework is certainly not useful for determining quality measures 
that are directly linked to the time evolution of the system, such as the question: "What 
is the probability that the network is connected for a time frame of length £ ?" I n order 
to answer such questions, we need to make specific assumptions on the motion of nodes. 

Certainly, it would be possible to incorporate one of the common explicit mobility 
models, like random waypoint or Brownian motion, into our context. From a more general 
point of view, however, these explicit models seem to be rather ad hoc and include some 
aspects that are not really motivated by properties of the real network system (such as 
discrete time steps). These technicalities could even complicate an explicit analysis more 
than necessary. Therefore, it might be desirable to consider a model with more natural 
assumptions, or explore and compare several such modelling alternatives. 

On the mathematical side, passing to such an analysis - without assuming a discrete 
time scale - would mean that we pass from our finite-dimensional sample space Q n to an 
infinite-dimensional space of functions. The base for such an ab origine calculation could be 
found in the theory of stochastic integrals and stochastic differential equations; while this 
field is well established |Pro03[ lPKL04j . it would certainly increase the technical complexity 
of our analysis by far, compared with the rather elementary mathematical methods used 
in the present work. Still, this might be a promising subject for future research. 



Appendix A 

Some Mathematical Machinery 



A.l The standard simplex in higher dimensions 

In our analysis, we often deal with a specific volume in R n , the n- dimensional standard 
simplex, defined as 

n 

K:={*G [0,1]" | ;>>,<!}. (A.l) 
i=i 

We are also often lead to the top surface of V n , which we denote by T n and define it as 

n 

T n := {x G [0, l] n | 5>* = 1}. (A.2) 

i=i 

T n is an (n — l)-dimensional manifold spanned by n corner points, which are all located at 
equal mutual distances; specifically, 7\ is a single point, T 2 a straight line, T 3 an equilateral 
triangle, and T 4 a tetrahedron. 

In this appendix, we discuss several properties of V n and T n , where it seems appropriate 
to develop them separately from the main text. Specifically, we calculate certain integrals 
over V n and T n that turn out to be important for our argumentation. 

First of all, let us introduce Heaviside's theta function 

6{x):={ 1 ifX ^°' (A.3) 
V ; [0 ifx<0; K ' 

note that 9(x) = 9(\x) for all A G 1R + , x G R. Also, we have 6(—x) = 1 — 9(x) except for 
x = (this set of zero volume can be neglected in integrals). Using the 9 function, we can 
express the integral of any function / over V n as follows: 

/ d n xf{x)= [ d n xf{x)9(l-^Txi). (A.4) 

JV„ i[0,l] n i= i 



We will now calculate certain integrals over V n explicitly. 
Lemma A.l. For any n G N and k G N , we have 

d n x (i-y Xl ) k = 7 J± 



I 

JVr, 
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Proof. We will prove the relation by induction on n. For n = 1, the proposition reads 

i 1 

dx (1 - xf = - , (A.5) 

v ' k + 1 v ' 

which is easily checked by direct calculation. Now let the proposition be true for n — 1 in 
place of n, with k e No being arbitrary. Using Eq. (|A.4I) . we calculate 



„ n „ „i n n—1 

/ dtx (l - ^ Xi) = / d n_1 a; / dx„,(l-^Xi) fl(l — — x ra ) 

= / d^xe^-J^Xi) *~ ^(l-Vii)'. (A.6) 

Setting a = 1 — X}i=i ^ the integral in x n can be elementary solved as 

f dxn (a - x n ) k = ^ [ - (a - x n ) k+1 ] a Q = -pja* 1 - (A.7) 
Inserting this result into Eq. ([A. 61) . we have 



n-1 

fc+1 



■>Vn j =1 -'Vn-l ' j = 1 

By induction hypothesis, this evaluates to 

(Jfe + 1)! fc! 



(A.f 



/ d n x (l-Yx t ) k = , 1 , {k 
Jv n K tt ] fc + l(n-l 



(A.9) 



+ Jfe + 1)! (n + Jfe)!' 

as desired. □ 

Our next task is to calculate similar integrals over the top surface T n of V n . For 
calculating such an integral of some function /, 



/ dS(x)f(x), (A.10) 



where dS(x) is the surface element of T n , we need a coordinatization of the surface and 
the length of its normal vector. Since the surface is characterized by the equation 



1-^^ = 0, (A.ir 



i=i 



coordinates are simply given by e.g. (x\, . . . , ar n _i) G V^_i, setting x n = 1 — YllZi x ii anc ^ 
the normal vector is easily seen to be (1, 1, . . . , 1) e M n , so its length is y/n. Thus, for 
n > 2, we can calculate the integral as 

„ „ n— 1 

/ dS(x) f(x) = Vn~ d n - 1 xf(x u ...,x n ^ 1 ,l-J2^i)- (A-12) 

JTn JVn-l i—i 
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For n=l, the surface T n is a single point, and we have f T dS(x)f(x) = /(l). We 
shall often represent the integral in a different way: Using the "delta valued measure" 
concentrated on T n , we can rewrite Eq. (jA.12|) as 

/ dS(x)f(x) = y/n f d n x5{l-y^x l )f{x). (A.13) 
Jr n J[o,i]» i=1 

which also holds for n = 1. In many situations, we prefer the latter form of notation, 
since it expresses the symmetry between the n different coordinates more directly. The 
reader unfamiliar with delta-valued measures |GS64] can always replace this expression 
with (jA.12j) if in doubt. 

We will now calculate some commonly used surface integrals. 

Lemma A. 2. For any n G N, k G No, and j G {1, . . . , n}, we have 

jo/ \ k \fnk\ 

dS{x)x k j = - — — 7 — . 

Tn 3 (n + k-l)\ 

Proof. The statement is easily checked for n = 1; so let n > 2 in the following. Due to 
symmetry reasons, we can choose j = n without loss of generality. Now setting f(x) = x\ 
in Eq. (JA.12j) . we see that 

„ „ n— 1 

/ dS(x)x k n = y/rl d n ~ x x (1-J2 x i) k - ( A - 14 ) 

J T n J Vn—1 i=l 

The integral on the right-hand side is known by Lemma fA.l| inserting that expression, we 
can immediately show the proposed result. □ 

Let us note some consequences of the previous lemmas: Setting k = in Lemma IA.11 
we can calculate the volume of V n as 

vol(K) = -.. (A.15) 

71 ! 

In the same way, setting k = in Lemma IA. 21 we can determine the (n — l)-dimensional 
volume of T„ as 



vol(T„) = (A.16) 

By the latter result, we can easily write down the probability measure of equal distribution 
on the surface T n , which fulfills dfi^~ eci (x) = (vol T n )~ 1 dS(x). We can summarize this as 
follows. 

Proposition A. 3. The measure of equal distribution over the surface T n has the form 

n 

u,l-^x) = (n-l)\5(l-Y. x *)> ( A - 17 ) 

i=i 

considered on the space [0, l] n . For any n G N, k G No, and j G {1, . . . , n}, we have 

-l 



62 



Appendix A. Some Mathematical Machinery 



The second part of the proposition follows directly from Lemma arid Eq. (|A.16J) . We 
now turn to another often-used relation, which might be described as a scaling argument 
on T n . To that end, let ey) = (0, . . . , 0, 1, 0, . . . , 0) (with the 1 in the j-ih place) denote 
the j-th standard unit vector in M. n . 

Lemma A. 4. Let n G N, j G {1, . . . , n}, and A G (0, 1), and let f : T n — ► R be integrable. 
Then 

[ d^(x) 6{x 3 - X)f(x) = (1 - A)"" 1 / dfi-*(x) /((l - \)x + Ae (i) ). 

J[0,l] n J[0,l] n 

Proof. Since the integration measure does not change when permuting the variables, we 
can assume without loss of generality that j = n. By Eq. flA.12|) . we have 

d^~ e \x) 9(x n - X)f(x) 

[0,1]" 

/n— 1 n— 1 
d n ~ l x 9(1 - Xi ~ A) f(xi, ar n _i, 1 - ^ x i) 

™-i i=i i=i 

„ n— 1 n— 1 n— 1 

= (n - 1)! / cf^x 0(1 -Vii) - A) - j><) /(rci, . . . ,x n -i, 1 -V4 
■W i= i i= i i= i 

(A. 18) 

Since (1 — A) < 1, the first theta function is redundant in view of the second one. Then, a 
variable transformation x\ = (1 — A) _1 Xj leads us to 



/ d^ e(x n - x)f( x ) 

J [o,i] n 

„ n— 1 

= (" - I)' (1 - A)"" 1 / ^"V /((l - A)xi, . . • , (1 - A)<_i, 1 - - A)^) 

= (1 - A)"- 1 / d^~ c \x') /((l - A)*' + Ae (n) ), (A.19) 

J[0,1]" 

which was to be shown. □ 

Using the previous lemma, we will establish a related technical result which turns out 
to be useful for our purposes. 

Lemma A. 5. Let n G N, n > 2, and let j, k G {1, . . . , n} with j ^ k; furthermore, let 
X G (0, 1). Then 

[ d/Z~ eq 0(X - Xj - x k ) = 1 - (1 - A) n - 2 (1 + (n - 2)A). 
J[o,i] n 
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Proof. In the case n = 2, both sides of the proposed relation vanish; so let n > 3. Without 
loss of generality, we may assume j — n — 1 and k = n. Observe that 

/ d^( x )e(\-x n ^-x n ) 

J [0,1] n 

/n— 1 
rf n - 1 x^(A-x n _ 1 -(l-^x,)) 
n-l j =1 

„ n-2 /-l-X]"-! 2 x i 

= {n-l)\ d n - 2 x6(\-l + J2 x i) *~ dx ^-i 

JVn-2 i= l JO 

„ n— 2 n— 2 

= (n - 1)! / d n ~ 2 x 9(X - (1 - x i)) (! ~ ( A ' 2 °) 
We can rewrite this expression as an integral over T n _\\ 

I dril~ eq (x) 6>(A - x n _i - x n ) = f drilzl q (x) 0(A - x n _i) x n _i 

7[o,i]" [ n ~ Z Y- ^[o,i]" " v ' 

l-0(a;n-l-A) 

= (n-l)( dfilzt q {x) x n _i - / du^ZT{x) #(x„_i - A) x n _i). (A.21) 

V 7[0,l]" ^[0,l] n / 

The first integral expression is known by Proposition IA.3| on the second one, we can apply 
Lemma IA.41 This yields 



f dri T n - e \x) 6{\ - x n ^ - x n ) 

J [0,11" 



I- in- i)(i _ A)"- 2 / df n ZT(x) (A + (1 - A)ar n _i). (A.22) 



Again applying Proposition IA.31 our result is 

dfil- e \x) 0(A - x n _! - x n ) = 1 - [n - 1)(1 - A)"- 2 (A + i— ^) 
[o,i]" n-l 

= l-(l-A)"- 2 (l + (n-2)A), (A.23) 
as proposed. □ 

A. 2 The inclusion-exclusion formula 

At several points in the main text, we make use of the well-known inclusion- exclusion 
formula, which allows us to calculate the probability of certain events easily. We formulate 
it here for reference. 

Theorem A. 6. Let Q be a sample space, and let Ai,...,A n C £1 be events in it. For 
k G N , let 

B k := {uj G VL I oj G Aj for exactly k values of j}, and 
Ck ■= {uj G Q | uj G Aj for at least k values of j}. 
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Then we have 

n 

P(B k ) = 

j=k 

n 

P(C k ) = X)( 

j=k 

where Sj is defined as 

Sj--= P(A mi n...nA mj ); 

{mi,...,mj } 

the sum runs over all subsets {mi, . . . , rrij} C {1, . . . ,n}. 

A proof of this formula can be found in most textbooks on elementary statistics - see, 
for example, the book by Krengel |Kre91| Sec. 3.4]. 



S 



3i 



A. 3 Statistical limits 

In this appendix, we will state some familiar limit theorems that are useful in our discussion; 
they are reproduced here for easier reference. The first of these is Stirling's formula, which 
gives an approximation of the factorial n\ for large n. Its precise form is: 

Theorem A. 7. For each n G N, there is a 6(n) G [ 12r [ +1 , such that 



ft! 



A proof can be found e.g. in (Kre91l Appendix to §5]. It follows in particular that 
ft! > (ft/e) n for all iieN; this is the inequality we will actually use. 

The next result which we want to note (in fact a consequence of Theorem I A. 7)1 is the 
Theorem of de Moivre-Laplace, which tells us about the convergence of binomial probability 
distributions to normal (Gaussian) distributions. To that end, let F be a random variable 
which is binomially distributed with parameters n and p; that is, 

Vz G {0, . . . , ft} : P(F = i) = ("Wl-p) 1- *. (A.24) 



We write a n = \/np{Y^j)). Further, let the function $ be defined as 



$(x) = — L / dy e- y2/2 ; (A.25) 



/2tt 

we additionally set $(— oo) = and $(+oo) = 1. Note that $(— x) = 1 — $(a;) for all x. 
The de-Moivre-Laplace theorem then states the following. 

Theorem A. 8. Let F be as above, and let a G R U {— oo} ; (3 G K U {+oo}, where a < (3. 
Then it holds that 

lim P(np + aa n < F < np + (3a n ) = $(/3) - $(«). 

n^oo 

For a proof, again see |Kre911 §5]. It is well known that the above theorem is only 
a special case of the more general central limit theorem; however, we shall only need the 
specialized form for our purposes. 



A. 4 A summation lemma 
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A. 4 A summation lemma 

This appendix presents an auxiliary result regarding a summation formula. The idea is to 
express the sum in question as a power series of a certain function, then using well-known 
relations for its derivatives in order to achieve the desired result. 



Lemma A. 9. Let j e N . Then 



k=0 



kiJ J 
k 



-1 ifj = h 
otherwise. 



Proof. For j = and j = 1, one checks by explicit calculation that the proposition is true. 
Now let j > 2, and let the function / be defined as 



/(*) 



:i-x) 



E 

k=0 



-X) 



Then we know from the expression on the right-hand side that 

j 



dx 



On the other hand, 



k=0 

dx 



k(3 7 )kx k - 1 
k 



-x 



) k k' J 



k=0 



x=l 



(A.26) 



(A.27) 



(A.28) 



since j > 1. Combining Eqs. (jA.27|) and (jA.28|) . we have proved the proposed result. □ 
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Notes on a Series of Publications by 
P. Santi et al. 



In a recent series of publications, P. Santi and D. M. Blough SI 502 S BQ3 . , as well as 



the same authors and F. Vainstein [SBV01J, have analysed the connectedness problem of 
MANETs using statistical models. Among others, they considered the very same mathe- 
matical model for a 1-dimensional MANET (with disconnected boundary conditions) that 
we have used in Sec. 13.31 The authors proposed a number of asymptotic estimates for the 
probability of connectedness; however, as mentioned in Sec. I3.3.2| these estimates are in the 
general case incompatible with the results of our analysis. The present author claims that 
several theorems established in [SBV01, SB02, SBQ3] do in fact not hold in the form stated 
there; this appendix will discuss counterexamples to those theorems, as well as pointing 
out inconsistencies in their corresponding proofs. 



In the following, we shall stick to the notation used in [SB03 rather than that used in 
the main text. This means in particular that we regard r and £ as two explicit parameters 
(rather than using the normalized radio range), that we will consider r and n as functions 
of £, and describe the limit of large systems as £ — > oo. (This is only a question of 
nomenclature.) 

Let us start with the upper bounds on the probability of connectedness as proposed in 
[SB03, Theorem 4]. The authors state the following. 

"Assume that n nodes, each with transmitting range r, are distributed uni- 
formly and independently at random in R = [0,£] and assume that rn = k£ln£ 
for some constant k > 0. Further, assume that r = r(£) <C £ and n = n(£) ^> 1. 
If k > 2, or k = 2 and r = r{£) > 1, then lim^^ P(CONN e ) = 1." 

(Here CONN^ is the event CONN-DB in our notation, and r -C £ means r/t — > 0, etc.) 
This statement is in conflict with our results: As a counterexample, consider r = £~ k In £, 
n = k£ k+1 , where k > 2. Then all conditions of the above statement are fulfilled; however, 
since Inn > (k + 1) ln£, one has 

k 

n r < £ Inn, (B.l) 

~ k + l v ' 

and thus -Pconn-db - > according to Corollary 13.121 
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In fact, the proof of Theorem 4 in the Appendix of [SB03 is inconclusive: After Eq. (2), 
the authors calculate the intermediate result 

2/ k In/ 2 
\nEMn, C)] < In - - — = In (B.2) 

where C = 2£/r, and fi(n,C) is a random variable whose details are not relevant here. 
Then they state: 

"If k > 2, or if k = 2 and r = r(£) ^> 1, then it is easily seen from this 
expression that lim nj c^ 00 lnE\p(n, C)] = — oo." 

However, this conclusion is not justified: In the above counterexample, one has 

2 2£ l+k ' 2 , x 

-> oo, (B.3) 



r ^/2-i inl 

thus it does not follow that the left-hand side of (jB.2|) converges to — oo. 

Note that the proof (and theorem) does hold in the case k = 2, due to the extra 
condition r ^> 1. It is also correct in the general case if one adds the condition that 
r > const, in the limit, or if one replaces the condition nr = k£\n£ with nr = k£\nn. (The 
proof can easily be adapted in the latter case.) 

The authors also presented a second, weaker result for the upper bounds SBV01, 



Theorem 4] , using a different proof technique. (The result is also reported within Theorem 3 



in |SB03 .) They claim the following: 

"Suppose n nodes are placed in [0, £} according to the uniform distribution. If 
rn G Q(£log£), then the r-homogeneous range assignment is a.a.s. connecting." 

(Here rn G Q(£\og£) means that £log£ = 0(rn), the r-homogeneous range assignment 
refers to the system considered above, and "a.a.s. connecting" means Pconn-db — > 1 in our 
notation.) This statement conflicts with our results as well, with a similar counterexample 
as above (where k is chosen sufficiently large). In fact, it is also in conflict with [SB03 ( 
Theorem 5]. The proof, as given by the authors, relies on Theorem 2 in |SBV01j . which 
reads: 



"Assume n nodes are displaced at random in [0,£\. Then, the probability that 
the r-homogeneous range assignment is connecting is at least 
1 - (£ - r)(l - f )V 

To see that this result is incorrect, remember that Pconn-db does not change when 
scaling both r and I together, i.e. when replacing £ with X£ and r with Ar, where A > 
is arbitrary. Exploiting this property, the above theorem leads to the conclusion that for 
any fixed £ and r, 

VA > : Pconn-db > 1 - X(£ - r)(l - - £ ) n ; (B.4) 
however, this would obviously result in Pconn-db = 1 for all parameter values. 
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The root cause of this error seems to be in the proof of the named theorem: Here, 
the authors define certain events DISCONNECTED^, where s G [0,£ — r] is a continuous 
parameter, such that 

DISCONNECTED^ (J DISCONNECTED^; (B.5) 

se[0,£-r] 

DISCONNECTED^ is the complement of our event CONN-DB. They then argue as follows. 

"An upper bound to P(DISCONNECTED£) can be derived by summing the 
probabilities P(DISCONNECTED^' r ) for all possible values of s. We thus have: 

pi— T 

P(DISCONNECTED^) < / P(DISCONNECTED*' r )ds [...]." 

Jo 

However, unlike the analogue case for a finite union of events, this is not a valid consequence 
of Eq. (jB.5|) - passing to the integral for the "summation of probabilities" is by no means 
justified. 

For the lower bounds on the probability of connectedness, Theorem 5 in jSB03j states: 

"Assume that n nodes, each with transmitting range r, are distributed uni- 
formly and independently at random in R = [0,£], and assume that rn = 
(1 — e)£\n£ for some < e < 1. If r = r(£) G Q(£ e ), then the communication 
graph is not connected w.h.p." 

Here "not connected w.h.p." corresponds to Pconn-db 1 in our notation. This 
theorem is compatible with the present work. In fact, using that the relation between ln£ 
and Inn is fixed by the requirement r G Q{£ e ), one can use Theorem 13.111 to show that 
-Pconn-db —* under the conditions given. 

Under more general conditions, Theorem 6 in jSB03j claims the following result: 

"Assume that n nodes, each with transmitting range r, are distributed uni- 
formly and independently at random in R = [0, £} and assume that r = r(£) £ 
and n = n(£) ^> 1. If rn £\n£, then the communication graph is not con- 
nected w.h.p." 

This statement again is incompatible with the results in Sec. 13.31 As a counterexample, 
let n = ln£ and r = £/lnln£, thus fulfilling all prerequisites of the theorem. In this case, 
we have 

£ \n£ \n£ 

nr = - — - — - = — — - — — £ \nn>2£ Inn for large £; (B.6) 
mhu (mku)^ 

so Corollary 13. 121 tells us that Pconn-db — ► 1- 

For its proof, the cited Theorem 6 of SB03J relies on |SB02[ Theorem 4] . The proof of 
that theorem, located in the Appendix of [SB02 , is in fact inconclusive: Defining C := £/r, 
the authors note 

"Observe that the condition ^ C ra C £\og£ implies that C < n < ClogC 
[...]." 
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However, in the general case, this implication does not hold: In the above counterex- 
ample, we have in fact 

t In I 

£<^rn = — — <£ln£, (B.7) 
In hxz 

but it follows from C = £/r = ln\n£ that 

n = ln£^lnln£ In In In ^ = C In C. (B.8) 



Thus, one cannot conclude n <ti C log C, and the subsequent arguments in |SB02] do not 
apply- 

In conclusion, let us briefly mention that Theorem 7 of |SB03j . which summarizes most 
of the propositions discussed above, does consequently not hold in the stated form. 



Index of Notation 



Asymptotic behaviour of functions. For two functions f,g, we write f = 0(g) or 
f(x) = 0(g(x)) if f(x) < g(x) ■ const in the limit being considered (usually x — > oo). 
The notation / = 0(g) is used as an abbreviation for / = 0(g) A g = 0(f). We write 
/ ~ g to denote that f(x)/g(x) — > 1. For sequences rather than functions, we use similar 
notation. The sign is used in a more qualitative sense in heuristic argumentation, 
meaning "approximately equal to" (in a sense to be specified later). 

Vector notation. Vectors (i.e. elements of some ¥L n ) are denoted by boldface symbols, 
while their components are denoted in normal typeface; e.g.: x = (x%, . . . ,x n ). We do not 
always explicitly specify the dimension of the underlying vector space where it is apparent 
from the context. 

Symbols and abbreviations. The following table lists symbols and abbreviations fre- 
quently used in the text, with a reference to their definition or first occurrence. 



symbol 



description 



reference 



CONN-DB 
CONN-PB 
k-DISCONN-DB 
k-DISCONN-PB 

E[F] 

M EV 

t 

n 

N 

No 

-Pev 

0(f) 



event of connected MANET with disconnected Eq. ()3.56|) 
boundary conditions 

event of connected MANET with periodic bound- Eq. (|3.2())1 
ary conditions 

event of fc-disconnected MANET with discon- Eq. (j3.21j) 
nected boundary conditions 

event of fc-disconnected MANET with periodic Eq. (|3.55|) 

boundary conditions 

j-th standard unit vector in W 1 

expectation value of a random variable F Eq. (12.5)1 

subset of Q n associated with an event EV Sec. 12.21 

spatial extent of the MANET Eq. (ETTJ) 

total number of network nodes in the MANET Sec. 12.11 

= {1,2,3,...} 
= {0,1,2,...} 

probability of an event EV Eq. ()2.7)) 

see "asymptotic behaviour of functions" above 



(continued on next page) 
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Index of Notation 



symbol description reference 

r radio range of a MANET node Sec. 13.11 

R+ ={xeR\x>0} 

M+ = {x e R I x > 0} 

T n top surface of V n Eq. ()A.2j) 

V n n-dimensional standard simplex Eq. (jA.lj) 



5( • )c/ n x delta-valued integration measure Eq. (|A.13j) 

Xev characteristic function of an event EV Eq. (J2.6|) 

1] limit of np — In n as n — > oo Thm. I3.8I 

Pn~ eq measure of equal distribution on T n Prop. IA.3I 

limit of np as n — >• oo Def. I4.ll 

p = r/£, normalized radio range Sec. I3.ll 

9( ■ ) Heaviside's theta function Eq. (jA.3|) 

0(/) see "asymptotic behaviour of functions" above 

fi n sample space for a MANET with n nodes Eq. (|2.1j) 

F = E[F], expectation value of a random variable F Eq. ()2.5|) 

[x] = max{/c G Z | /c < x}, Gauss bracket of x 

M c complement of a set M 
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